A FEW COMMENTS ON CLASSICAL ELECTRODYNAMICS 



K. Ghosh* 

Department of Engineering Science, Haldia Institute of Technology, 
Hatiberia, Purba Midnapore - 721 657, India. 



I. ABSTRACT 

In this article we discuss a few aspects of the space-time description of fields and particles. In section 
II wc will show that a line is not just a collection of points and we will have to introduce one-dimensional 
line-intervals as fundamental geometric elements. Similar discussions are valid for area and volume-elements. 
In section III we will discuss the fundamental nature of fields. In section IV we will discuss non-equivalence 
of the Schwarzschild coordinates and the Kruskal-Szekeres coordinates. In App.B and App.C we make a 
comparative study of Quantum Field Theory and Quantum Mechanics and disagreements associated with 
probabilistic interpretation of these theories with space-time dimensional analysis. In App.D and App.E we 
discuss the geometry of Robertson- Walker model and electrostatic behavior of dielectrics respectively. In 
Supp.I we discuss the regularity of Spin-Spherical harmonics. In Supp.II we will give a causal description 
of interference and diffraction in terms of phase velocity of waves. In Supp.III we will discuss the regularity 
conditions required for the integral version of Gauss's divergence theorem and rederive Gauss's law in Elec- 
trostatics. We also show In Supp.IV we derive the complete equation for viscous compressible fluids and 
make a few comments regarding some contradictions associated with boundary conditions for fluid dynamics. 
In Supp.V We make with a few questions in We make a few comments on seeing by human brains in Supp.VI. 
We make a few discussions on radiation in isothermal enclosures and uncertainity principle in Sup. VII. In 
Supp.VIII We make a few discussions on Newton's laws for conservative force fields. In Supp.IX we derive 
an expression for scalar field entropy in the Schwarzschild black-hole background using the brick-wall model 
which is different from that of ' Hooft. We conclude this article with a few discussions on surface integral of 
vector fields in the last section " Miscellaneous" . 
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II. GEOMETRIC CONTINUITY AND COORDINATIZATION 

In this section we will review a few basic concepts of geometry and coordinatization. We will first give 
a precise definition of geometric continuity based on the concepts of homogeneity and extension. We will 
demonstrate that a line-element, a geometric continuum of one-dimension can not be a collection of points. 
Line-elements, area-elements and volume-elements are as fundamental as points. It follows naturally that 
an area-element can not be described as a collection of line-elements. This leads to a reformulation of the 
definition of dimension for a geometric continuum. We will consider the fact that real numbers form a dense 
set: for real numbers only the concept of "two numbers separated by an interval" or "two numbers are same" 
have meaning but "two numbers adjacent" is not defined. We will illustrate the relationship of these two 
aspects in the context of coordinatization. We will also state a corresponding version of the completeness of 
the real numbers. The definition of geometric continuity as given in this article raises the question whether 
time is a geometric continuum. We will briefiy address this issue through the conceptual foundations of the 
General Theory of Relativity. We will also discuss in brief the the relationships of these aspects and the 
three spatial dimension with the the kinematical and dynamical properties of the Elementary particles and 
Electromagnetic fields. We will conclude this section with a few consequences. 

In this section we will demonstrate that a line-element is not a collection of points. Line-elements, area- 
elements and volume-elements are as fundamental as points. We will consider the proposition similar to the 
fact that rational numbers form a dense set: for real numbers only the concept of "two numbers separated 
by an interval" or "two numbers are same" have meaning but "two numbers adjacent" is not defined. We 
will illustrate the relationship of these two aspects in the context of coordinatization. 

The concept of geometric continuity should be based on two aspects: extension and homogeneity (more 
properly uniformity). An elementary particle, a tree if we neglect microscopic non-homogeneities, space and 
space-time (we will illustrate this example later) are examples of geometric continuum. In the most broad 
sense geometric continuity is distinguished according to dimension. An elementary particle, an electron, is a 
three dimensional geometrical contimium. If we consider a plane passing through the center of the electron 
(an equatorial cross-section) and shrink the transverse geometric extensions to zero we get a two-dimensional 
geometric continuum: a disk. If we consider a diameter of the disk and shrink the transverse extensions 
along the disk to zero we get a one-dimensional geometric continuity: a line. This definition of dimension 
is different from the conventional definition in terms of possible linear motion of point particles [32] but is 
proper as we will illustrate later. A point is a geometric object with zero extension. Thus whatever may be 
the amount of content (cardinality), a collection of points can not give a geometric continuum, a geometric 
object with finite homogeneous extension. We can illustrate this aspect in the context of the infinite straight 
line, a geometric continuum of one dimension. The one-dimensional homogeneous extension of a line-element 
on the straight line is characterized through the concept of length. Qualitatively, length characterizes finite 
homogeneous extension in one-dimension and also ordering of points with respect to a given point chosen 
as origin. Quantitatively, length gives a coordinatization scheme characterizing line-intervals with a given 
line-interval chosen as unity. We will illustrate this aspect later. Thus, whatever may be the cardinality, 
a collection of zero-length points can not give a finite-length one-dimensional line. Also a homogeneous 
collection of points give us a single point as the length of a point is zero. We can illustrate this aspect in the 
following way: 

We consider a homogeneous linear array of marbles touching each other. If we now shrink the volume of 
each marble keeping them in contact (so that the array is always homogeneous) then in the limit that the 
volume of the marbles is zero we will get a single point. 

These discussions lead us to conclude that a lin(;-elenic!nt is not just a colk;ction of points but is a fim- 
damental geometric object. Similar arguments lead us to conclude that area-elements (two-dimensional 
geometric continuums) and volume-elements (three-dimensional geometric continuum) are fundamental geo- 
metric objects and we can not obtain points without breaking the continuum geometry. To illustrate, let us 
consider a square, a geometric continuum of two dimension. We can not consider the square to be a collection 
of straight lines parallel to the base as each straight line has zero transverse extension and arguments similar 
to as above apply. 

We will now discuss the concept of coordinatization which is an ordering process through length with 
respect to a point chosen as the origin. In the conventional scheme the infinite straight line is considered 
to be a collection of points and real numbers should have a one-to-one relationship with the points to have 
a proper ordering. In the context of the discussions in the preceding paragraph this scheme can not work 
as we can not arrange, according to length from a given point, a homogeneoiis collection of points each of 
which is of zero length. We have also illustrated earlier that a homogeneous collection of points give us a 
single point. A plane is a two dimensional object and any point on the plane can be characterized, as we 
will illustrate, through the choice of two intersecting lines on the plane and constructing a coordinate system 
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defined in terms of intervals witli respect to tlic point of intersection. 

We now consider the concept of coordinatizationmorc precisely. As we have discussed earlier line intervals 
as fundamental geometrical entities to form lines which are geometric continuum of one dimension. 

We define a dense set to be a set with an ordering such that there exist an element between any two 
distinct elements of the set. We will now prove that the real numbers form a dense set [33,34]. That is for 
real numbers only the concept of "two numbers separated by an interval" or "two numbers coincident" have 
meaning but "two numbers adjacent" is not defined. That is for any two real numbers we can always find a 
real numbers lying between the two. 

Firstly, the two representations of 0.1, 0.0999 (where all the decimal places starting from the 2-nd are 

9) and 0.1 are same as any number added to 0.0999 gives a number greater than 0.1 and any number 

subtracted from 0.1 gives gives a number less than 0.0999 This two representations also gives us the 

same result if add/subtract a number with/from both the numbers. 

An arbitrary real rational number may be expressed as r = n.xiX2X3 Xp, where n is an integer, 

Xi = 0, ...,9 for i < p, Xj — for j > p, p may be arbitrarily large but finite and Xp ^ 0, i.e., the 
sequence of the decimal places is finite. Whereas an irrational number and a recurrent decimal is given 

by: ir = n.x\X2Xz where all the decimal places are significant. We will consider r,ir > in the 

following section. 

We now prove the following proposition: 

For real numbers only the concept of "two numbers separated by an interval" or "two numbers coincident" 
have meaning but "two numbers adjacent" is not defined. That is given any two real numbers we can always 
find a real numbers lying between the two. 

If r is an integer (> 0) and M = r.xiX2X3 Xp, with all but Xp (= 1) are zero and p may be arbitrarily 

large, then the number of decimal places having value zero before Xp = I can be increased arbitrarily to 
construct a number less than M and closer to r and it is not possible to define a number to be adjacent to r. 

Similar will be the case for N = {r — 1).999 Xi ,with all the decimal places up to Xi are 9 for i arbitrarily 

large, to construct a number greater than N and closer to r, i.e, the number of decimal places having value 
9 can be increased arbitrarily and it is not possible to define a number to be adjacent to r. We can illustrate 
this issue with the representations discussed before the proposition. The numbers r and (r— 1).999... where 

all the decimal places are 9 equivalent. Thus we can form the number (r — 1).999 Xi = 899.. and applying 

arguments similar as above to the i-th decimal place we can show that it is not possible to define a number 

adjacent to r. This is also same as subtracting M ~ 0.000 a;=l from r and applying arguments similar 

to those discussed at the beginning of this paragraph. 

Similar as above will be the arguments with (n > 0) r = n.xix^xz Xp, M = n.xiX2X3 xj (a;, 

same for both for 1 < i < p) with all Xi {i > p) but Xj (= 1) (j > i) are zero for j arbitrarily large and for 

N = n.xiX2X3 j/p999....Xj with the p -th decimal place yp = Xp — 1 and all Xj for j arbitrarily large 

are equal to 9. 

For irrational numbers like (ir)i = \/2, tt we can prove the proposition using the rational approximation 
as considered in the standard literature. We can also apply the above arguments for recurrent decimals. 
Following the constructions in the preceding paragraphs it is easy to show that between any two irrational 
number there exist a rational or irrational number. 

For r < the corresponding arguments to prove that "two numbers are adjacent" is not defined are very 
similar as in the preceding paragraphs. 

These arguments together with the completeness of the Decimal number system (in the context of number 
systems) prove the above proposition. 

We now consider a few aspects about irrational numbers defined in the conventional way. If we represent 
irrational numbers in terms of nested families of closed intervals [32] with rational boundaries then we will 
have to assign a corresponding non-zero length (as for Euclidean space length is same as the difference of the 
rational boundaries) to the null rational interval representing a single point on the real line corresponding 
to the irrational number and this is not well-defined. On the other hand we can assume that between two 
rationals there exist a set of irrational numbers with forming a line- interval with continuum cardinality [6] . 
(This has some consequence for first countability of the real line as we will discuss at the end of this section) 
We now consider the null rational interval [0, r]. Whatever r can be we can multiply this interval with a 
factor of the form (10)*^ (r can be a surreal number and M can be a supcrreal number [35]). The resultant 
interval should contain an infinite number of rationals [33] although all the numbers in the interval [0,r] 
(apart from the two boundaries) are irrational and multiplying an irrational number by (10)^ can not give 
us a rational number. Thus we take the point of view of [36], i.e, the rational approximations of an irrational 
number form a dense set about the irrational number and the rational approximation give us the irrational 
number in the limit when the number of significant decimal places are infinite, i.e, when the number of 
significant decimal places never terminate. For a recurrent decimal and an irrational to be adjacent the 
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irrational should differ from the rational at some finite decimal place and arguments similar as discussed 
earlier demonstrate that it is not possible to define an irrational and a recurrent decimal to be adjacent. 

These discussions also lead us to conclude that it is not possible to define a successor for real numbers. If 
we do not accept the point of view of [34] to break the fact that real numbers form a dense set then this aspect 
of real numbers is in accordance with geometric continuity as defined at the beginning of this article. That is, 
as far as coordinatization is concerned, the real numbers characterize line-intervals, which can be; arbitrary, 
with respect to a point chosen to represent the value zero, the origin. This scheme of coordinatization diflFers 
from the conventional way where we consider the the infinite straight line to be a collection of points with 
continuum cardinality and try to have a one-to-one relationship between the points and the real numbers. 
In this context we can reformulate the completeness of the real numbers in the following way: 

If we chose any interval on the infinite straight line (the Real line in the sense that the intervals on the line 
is characterized by real numbers with a point on the line chosen to represent the number zero) the length 
can be expressed either in terms of a rational number or an irrational number. Similar will be the situations 
with the higher dimensional volume elements. 

We note as the real line cannot be defined as a collection of points we have to introduce the concept of 
collection of one-dimensional line elements as a fundamental mathematical entity which can be characterized 
into four classes: 

(i) one-dimensional line elements without boundaries 

(ii) one-dimensional line elements with one boundaries 

(iii) one-dimensional line elements with two boundaries 

(iv) one-dimensional line elements with the two boundaries identified (e.g, a circle) 

We now demonstrate the consistencies of the discussions regarding geometric continuity at the beginning 

of this article in the following way. We can define the radius of a circle as the perimeter over 2tt. As is 
proved earlier the values of the radius can only be defined through intervals and the concept of two concentric 
circles with adjacent values of radii is not defined. This feature is consistent with the fact that we can not 
obtain a two-dimensional disk from a collection of one-dimensional circles. We have to introduce the two 
dimensional circular strips togetrher with the point origin as fundamental mathematical objects to construct 
a two-dinic!nsional disk in the above way. Similar arguments for two dimensional spheres (where the radius 
is now defined as the positive square root of the area of the two-dimensional sphere over 47r) illustrate that 
we can not have a three dimensional volume element from a collection of two dimensional spheres and we 
have to consider three dimensional volume elements as fundamental mathematical entities together with the 
point origin. These discussions can be extended to higher dimensions. 

A natural consequence of the above discussions is the fact that even if a physical variable is a geometric 
continuum its values are defined through intervals with respect to a reference value and form a dense set. 
Existence of quanta are restrictions on the values of these intervals leading to isolated speetriim. 

The finite volume of the elementary particles and their displacements are in accordance with the geometric 
continuum structure of space. Within the context of classical mechanics time is realized through the changes 
in the configuration of a system. To illustrate we can consider the coordinates of a moving particle to 
characterize time. In the conventional definition of a line as collection of points we can define time to 
be a geometric continuum, as a collection of points, by a one-to-one correspondence between the possible 
coordinates of the particle on the line and the instants of time. However as we have discussed in this article a 
line is not a collection of points. Line-intervals are fundamental geometric objects. Only changes of position 
of the particle by a non-zero line-element is meaningful. The coordinates do not form a geometric continuum 
but form a dense set characterizing line-intervals with respect to a chosen point as the origin. Thus if we 
define the instants of time through a one-to-one correspondence with the possible coordinates of the moving 
particle we will have a dense set structure for time. However the General Theory of Relativity is a theory 
governing the dynamics space-time intervals themselves. The kinematical equivalence of space and time 
intervals indicates that time should be a geometric continuum in the same way that space is. The space-time 
coordinates of events form a dense set. 

Unlike the motion of points, line-elements and area-elements the motion of three-dimensional objects 
form a geometric continuum of the same dimension as that of the objects themselves. Thus the three 
dimensional spatial geometry of the Universe is realized through the finite three dimensional volume of the 
fundamental particles and the finite three dimensional volumes of the fundamental particles lead to three 
spatial dimensions for the Universe. Dynamically, there are two possible motions on a plane, translations and 
rotations. Dynamically infinitesimal rotations can be uniquely characterized through their magnitude and 
the direction of the normal to the plane of rotation chosen through proper convention. In three dimensions 
the normal to the plane of infinitesimal rotation is uniquely characterized up to orientation, determined 
by proper convention, and is contained within the three dimensional spatial geometry. Higher dimensions 
greater than three break the uniqueness of the normal to the plane of rotation. Thus three dimensional spatial 
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geometry is self -complete to describe the dynamics of matter particles uniquely. The quantum mechanical 
description of spin and polarization of the elementary particles is also significant in this respect. Also the 
concept of orientation with a proper convention is an essential aspect to formulate the laws of Classical 
Electrodynamics consistent with the energy conservation law [24]. The fact that in general, apart from the 
case when the normal to the plane of the loop is perpendicular to the magnetic field, a time- varying magnetic 
field always induces an electric current in a closed loop obeying Faraday's law (consistent with the energy 
conservation law) indicates that spacial geometry is of three dimensions as far as Classical Electrodynamics 
is concerned. 

To summarize, in this section we have established that a line-interval is not a collection of points and 
we will have to consider line-intervals as fundamental geometric objects. Similar discussions are valid for 
higher dimensional geometric continuum. We have used the fact that real numbers form a dense set: "For 
real numbers only the concept of "two numbers separated by an interval" or "two numbers coincident" 
have meaning but "two numbers adjacent" is not defined" and discussed the relation of these two facts in 
the context of how to define coordinatization. We have also discussed in brief the the relationships of these 
aspects and the three spatial dimension with the the kinematical and dynamical properties of the Elementary 
particles and Electromagnetic fields. 

We conclude this section with a few discussions: 

Firstly, we can not express the Cartezian plane in terms of the direct product R^xR^. R^ is not just a 
collection of points. Hence it is not possible to express R^ as the union of the corresponding number of R^s. 
Also we have discussed earlier that we can not express a two-dimensional geometric continuum as a homo- 
geneous collection of one-dimensional geometric continuum. Since R-^ can not be expressed as a collection 
of curves, the tangent vector space at any point on R^ can not be identified with i?^. The above discussions 
are also valid in higher dimensions. This aspect also illustrate that, we can not define higher dimensions 
greater than one in terms of possible motion of a point particle as neither a two-dimensional area-element 
nor a three-dimensional volume-element is a collection of one-dimensional line-elements. These discussions 
are significant for Difi'erential Geometry. We now illustrate in brief the significance of these discussions in 
the context of classical statistical mechanics of a system of particles. We define the configuration-space to be 
the collection of the spatial coordinates of the particle. As we have illustrated earlier the configuration-space 
forms a dense set (a dense vector space when the dimension is greater than one). The momentum of a particle 
moving along a line is defined through associating possible values of the momentum with the possible posi- 
tions of the particle in a given time starting from a fixed point. For a single particle moving in one-dimension 
the corresponding momentum at any point is not a geometric continuum but form a dense set as changes in 
the position of the particle through arbitrary but non-zero intervals are only meaningful. This discussions 
are similar to those made for time. Also the discussions at the beginning of this paragraph indicates that 
in higher dimensions the momentum of the particle at any point, tangent vector to the possible trajectories 
of the particle at that point, do not form a geometric continuum but is a dense vector space. Thus the 
dynamical states of the particle is characterized by the coordinate dense set (the configuration-space) and a 
dense collection of the momentum space dense set. The fact that the coordinates and the momentums form 
a dense (vector) space indicates that physically we can only approximate the corresponding phase-space to 
a geometric continuum. The above discussions can be extended when the number of particles is greater 
than one. We will later illustrate the significance of these discussions in the context of the corresponding 
path-integral formulation of statistical mechanics. In passing we note the following: if we express a geometric 
continuum as a collection of points the cardinality and hence the amount of accessible microstates remains 
to the same irrespective of dimension. It is through coarse graining that one obtains different expressions 
for entropy depending on dimension. 

Secondly, let us consider two sets each containing a single object: a closed line-element. The line-elements 
are intersecting but non-coincident everywhere. The intersection of these two sets is a set containing points 
which are fundamentally different from line elements. Also as a line-element is not just a collection of 
points a set-theoretic union to form a closed line-element out of a set of half-open line-elements and a set of 
points is not well-defined as the above mentioned sets contains elements that are geometrically and thereby 
intrinsically different. Similar discussions remain valid in two dimensions when two area elements intersect 
along a line-element and in three dimensions when two closed volume-element intersect along a common 
boundary. These problems may be solved by defining a universal set containing all possible geometric- 
elements including points. The union or intersection of any two elements from this universal set give us 
elements which belong to this universal set. 

Lastly we discuss a consequence of the discussions on irrational numbers as discussed earlier. 

To illustrate we make some comments regarding the first and second countability of R^ discussed in 
Appendix: A of Wald [3]. A topological space X is first-countable if for every point p belonging to X there 
is a countable collection of open sets such that every open neighbourhood of p contains at least one member 
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of the collection. Whereas X is second countable if there is a countable collection of open sets such that 
every open set can be expressed as an union of open sets from this family. For R^, the open balls with 
rational radii centered on points with rational coordinates form such a countable collection open sets. When 
defined as in [6] there are infinite number of irrational numbers between two rational numbers. As far as 
first-countability is concerned, an open set V centered on a rational number with deleted peremeter on one 
of these neighbouring irrational number cannot contain opcin balls with rational radii centered on points 
with rational coordinates. Similar discussions can also be extended to i?". 



III. FUNDAMENTAL NATURE OF FIELDS 

Gravitational interaction which is universally attractive is described by the General Theory of Relativity, 
a theory of fields realized through the existence of space-time and some of its geometric properties, e.g, the 
curvature. The principle of equivalence leading to the fact that the guinea and the feather fall the same way 
in vacuum together with the fact that gravity violets the elementary quantum mechanical principles (as will 
be evident from this article) indicates that fields are as fundamental as particles. Fields and particles with 
their common and contradictory kinematical and dynamical features are the two fundamental constituents 
of nature. The attraction of opposite charges can only be explained through accepting the electric field 
as a fundamental entity of nature and not through only interaction mediating particle interaction (be the 
interaction mediating particles interacting with the sources or among thcimselvcs) although the repulsion 
between like charges may be explained in the later way with photons as the interaction mediator. Quantum 
mechanics is an incomplete approach to explain the Solar system microscopic physics purely in terms of 
wave-like properties whereas Quantum Field Theory is an incomplete approach to explain the Solar system 
microscopic physics purely in terms of particle-like properties. For the vacuum polarization explanation 
of the potential [18] in the Bhaba scattering process it is not obvious how, in terms of photon exchanges, 
the loops in vacuum will screen the charges of e~,e~^. Also an electron-positron loop with two external 
photon lines can not explain charge screening in any process as the particle-antiparticle pairs are created 
and annihilated after the interaction mediating photon had been created at one real particle and before 
the interaction mediating photon has interacted with the other real particle. These two theories are also 
contradictory to gravity (a theory of space-time) in many aspects as will be manifested in this article. 

During the last few decades a lot of efforts had been devoted to unify the general theory of relativity 
(describing the gravitational interaction) with quantum mechanics (describing the microscopic interactions 
of the elementary particles). Yet the conventional theory of quantum mechanics, based on unitarity and 
symmetries, is contradictory with general relativity in many respects, e.g. , formally infinite zero point 
energy associated with canonical quantization scheme, ultraviolet divergent energy density associated with 
vacuum fluctuations [1] for collapsing physical systems, unitarity violation for black hole evolution. 

To unify these two descriptions of nature we can proceed along two directions [2] : 

For the Euclede school space-time geometry is an abstract concept which exist irrespective of matter fields. 
For the non-Euclede school space-time does not exist independent of matter fields. Space-time is form of 
existence of matter and can not be conceived without matter. This feature is even more transparent from 
the facts, among others, that the universe is compact and there is no well-defined stress-tensor for gravity 
which interwinds matter fields and space-time geometry [3]. 

In the context of the general theory of relativity the conflict between the two schools arise in the following 
way: 

As soon as one derives the geodesic deviation equation from the principle of equivalence one "can" forget 
the source, ascribing the relative acceleration of the nearby geodesies to the space-time manifold. This, in 
contradiction to the philosophy of general relativity, may lead to think that the space-time manifold is more 
fundamental leading to the concept of quantum gravity irrespective of existence of the corresponding sources 
that will produce the fluctuating geometries (no source indicates no space-time geometry and observationally, 
quantum fluctuations in matter fields are negligible to produce significant alternations of the space-time 
geometry). 

Quantum gravity also led to many space-time geometries which are physically non-existent. One such 
example is the extreme Reissner-Nordstrom black hole which cannot be obtained through any realistic 
gravitational collapse [ ] . It also has vanishing Hawking temperature. This may be interpreted through the 
fact that Hawking radiation through pair production near the black hole event horizon is not possible as the 
metric do not change its signature across the horizon. 

We should keep in mind that the Einsteins equations in the regime of its validity determine the space-time 
geometry: the geometry of space-time as a whole is determined by corresponding matter fields described 
either in terms of some classical models or by a proper quantum theory. This gives a particular cosmology 
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(the closed or mathematically more properly compact universe picture) and if we view cosmology as a whole 
there is really no test body. 

We will now consider some aspects of the black hole spacc-tirnc geometry. In the process of gravitational 
collapse an event horizon, the black hole event horizon, is formed breaking the global CP invariance and 
giving rise to the Kerr- Newman families of black holes (the no hair theorems). The black hole event horizon 
may be defined as the causal boundary of the set of complete time-like geodesies which originates at the 
past time-like infinity and terminate at the future time-like infinity as classically nothing can come out off 
the horizon. The black hole space-time is usually described either in terms of the Schwarzschild coordinate 
system or in terms of the Kruskal-Szeckers coordinate system. In the Schwarzschild coordinate system the 
black hole event horizon is a two dimensional fixed point set of the time-like Killing vector field across which 
some of the metric components change sign. In the Kruskal-Szeckers coordinate system the event horizon is 
a two dimensional null surface across which the square of some of the coordinates change sign. We will now 
consider the non-equivalence of the two coordinate systems in detail. 

IV. NON-EQUIVALENCE OF THE SCHWARZSCHILD AND THE KRUSKAL-SZEKERS 

COORDINATE SYSTEM 

The Schwarzschild space-time is a Lorentz signature, static spherically symmetric solution of the Einstein 
equations when the Ricci tensor vanishes. This sohition describes the exterior geometry of a static spherically 
symmetric star and has been used to verify the predictions of general relativity for the Solar system. 

A space-time is said to be static if there exits a space-like hypersurface which is orthogonal to the orbits 
of the time-like Killing vector field. A space-time is said to be spherically symmetric if the space-like 
hypersurfaces contains 5*0(3) as a subgroup of the group of isometrics. The orbit spheres of SO{3) are 
isometric to the unit two sphere. These features together with the condition of the asymptotic Newtonian 
limit give the well-known Schwarzschild solution in the spherical polar coordinates [3]: 

ds^ = -(1 - 2M/r)df + {1- 2M/r)-^dr^ + r^[de^ + sin^ edO"^] (1) 

According to the Birkhoff's theorem [19] all spherically symmetric solutions with Rab = are static and 
the Schwarzschild space-time is the unique static spherically symmetric solution, up to diffeomorphisims, of 
the Einstein equations with Rab = 0. 

The norm of the time- like Killing vector field and (Vr)° in the orthonormal coordinates vanishes and some 
of the metric components are not well-behaved at r = 2M in the Schwarzschild coordinates. The proper 
acceleration of the constant r observers can be obtained from the geodesic equations in the Schwarzschild 
coordinates. This acceleration, a = (1 — 2M/r)~^/^M/r^, is divergent at the horizon (r = 2M). 

The ill-behavednes of the Schwarzschild coordinates is not a coordinate singularity like that of the spherical 
polar coordinate system where the azimuthal angular coordinate become ambiguous at the poles. All the 
ill-behavednes of the Schwarzschild coordinates at the horizon originate from that of the space-time metric. 
The curvature scalars calculated from the metric are well-behaved at the horizon unlike r = where the 
curvature scalars diverge. For ordinary stars this metric singularity at r = 2M is irrelevant as it is inside the 
star and the Schwarzschild solution is not valid in the matter filled interiors. However it is well-known that 
sufficiently massive stars can undergo gravitational collapse to form black holes and the metric singularity at 
the horizon is important. Several coordinate systems had been introduced to remove the metric singularity 
and to extend the Schwarzschild space-time where the Schwarzschild coordinate system is referred to covering 
a proper submanifold of the extended space-time. The metric in these c;xtendc;d c;oordinate systems are well- 
defined every where apart from the space-time singularity. The most well-known extension is the Kruskal- 
Szekers coordinate system. In this article we perform a comparative study of these two coordinate systems 
and show that they are not diffeomorphically equivalent. 

In this section we will follow the abstract index convention of Wald [3] and extend its significance in 
Appendix: A. 

According to the theory of relativity \i (p : M ^ M \s diffeomorphism then {M,gab) and {M,(l)*gab) 
represent the same physical space-time. Let a coordinate system x'^ cover a neighborhood J7 of a point p 
and a coordinate system y'^ cover a neighborhood V of the point (j){p)- Now we may use (p to define a new 
coordinate system x'^ in a neighborhood O = by setting x'^ = y^[(l>{q)] for q belonging to O. We 

may then take the point of view as </> leaving p and all tensors at p unchanged but inducing the coordinate 
transformation x^^ x"' . For to be a diffeomorphism should be non-singular [3,15]. According to this 
point of view two coordinate system covering a space-time can be taken to be equivalent if the corresponding 
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transformation coefficients arc not singular in tlieir common domain of definition otherwise an arbitrary 
smooth, function defined in one coordinate system may not remain smooth in the other coordinate system. 
To extend the Schwarzschild coordinate system one considers the two dimensional r — t part: 

= -(1 - 2M/r)dt^ + (1 - 2M/r)-'^dr^ (2) 
The Regge- Wheeler coordinate system is defined through the nuU-geodesics and is given by: 

n=r + 2Mln{r/2M - 1) (3) 
in this coordinate r — )• 2M corresponds to r* — )• — oo. The null coordinates are defined as: 

u = i - r* , = t + r* (4) 
A regular metric is obtained through the following transformation, 

;7 = -e-"/^^, ^ = 6"/^^ (5) 
The metric in these coordinates becomes: 

= ^udV (6) 

r 

As there is no longer a coordinate singularity at r = 2M (i.e at U = or V = 0) one extends the 
Schwarzschild solution by allowing U, V to take all possible values. However the transformation coefficients 

dU/dr = -d[{r/2M - iy/^e-^]/dr and dV/dr = d[{r/2M - Xfl'^e^ydr are singular at r = 2M 
and the extension is not diffeomorphically equivalent. Consequently as discussed at the beginning of this 
section the Schwarzschild coordinate system and the ([/, V) coordinate system do not represent physically 
the same space-time manifold. Consequently, according to Birkoff's theorem, the space-time represented 
by the {U,V,9,(j)) coordinate system is not a solution of the Einstein equations for a spherically symmetric 
black hole. 

Similar discussions arc valid for the Kruskal-Szekers coordinate transformations which are obtained 
through the following transformations: 

T = {U + V)/2, X = iV-U)/2 (7) 

and the metric becomes, 

ds^ = ^^^^ {-dT^ + dX^) + r\de^ + sin^ edcf") (8) 

The relation between the (T, X) and the (t, r) coordinates are well known and in the physical regions of 
interests are given by [4], 

X = (r/2M - lf/^e'''^'^cosh{t/4.M) (9) 
T = {r/2M - l)^/2e''/^^sinh(i/4M) (10) 

valid for r > 2M, and 



T = (1 - r/2M)V2e'-/4^cosh(i/4M) 



(11) 
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X = (1 - r/2M)i/2gr/4Mgjj^j^^^/4M) (12) 

valid for r < 2M. 

Again the transformation coefficients are not defined on tlic horizon and the Kruskal-Szekers coordinates 
do not give a proper diffeoniorphic extension of the Schwarzschild coordinate system. Hence the Kruskal- 
Szekeres coordinates is not a solution of the Einsteins equations for a spherically symmetric black hole 

The Kruskal-Szekers coordinate system had been introduced to eliminate a particular singular function (the 
metric components) in the Schwarzschild coordinate system through a singular coordinate transformation. 
This does not ensure that all singular tensors can be made regular in the new coordinate system and also 
tensors which are regular in the {t, r) coordinates can become singular in the (T, R) coordinates. To illustrate 
these features we consider the implicit relations between the two coordinate systems [1]: 

(r/2M - l)e''/^^ = - (13) 



2M = '<X^^ ^''^ 
The horizon in this coordinates are defined as X = ±T. 

Firstly the proper acceleration of the curves in Kruskal-Szecker's coordinate system which correspond to 
the constant r observers in the Schwarzschild coordinate system is given by a = {X'^ — T^)^^l'^\eJ' 1"^^ M lr^\. 
This is also divergent on the horizon. 

Secondly we consider the vector (^)", i? the proper rate of change of the curvature scalar R obtained 
from {dK)"" and the proper distance ds [i.e. the vector {^)[-§^\ Appendix:A]. The norm of this vector in 
the Schwarzschild coordinate system is (clRj dr)"^ and is finite on the horizon. Whereas the corresponding 
quantity in the {T^X) coordinates can be obtained from the following relations [apart from normalizing 

factors: (^), (f^) = [^]^/^ : 

dR _ dR dr dR _ OR dr 

dX~'d^dX' dT~'d^df 

and from equ.(13), 

dr SM^Xe-'-/^^ dr SM^Te-'^/^M 



dX r ' dT 



(16) 



and we have \R°-ks? = e^M^^^ [^f[x'^ - T^] = on the horizon although the r -dependent multi- 
plying factor in front of the Kruskal-Szecker's metric is finite at r = 2M. 

The unit space-like normal vector to the r = constant surfaces, which can be defined apart from r = 0, 
fc" = (^)" has unit norm (fc°fca = 1) on r = 2M although A;" — >■ as r — )• 2M which for an outside observer 
(r > 2M) may be interpreted as nothing can propagate radially outward at r = 2M. 

For two metric spaces the definitions of continuity is as follows [16]: 

Let {S, ds) and (T, dx) be metric spaces and let / : S ^ T be a function from S to T. The function / is 
said to be continuous at a point p in 5 if for every infinitesimal e > there is an infinitesimal 5 > Q such 
that 

dxlf {x) , f {p)] < e, whenever ds[x,p] < 6. (17) 
If / is continuous at every point of S then / is continuous on S. 

The definition is in accordance with the intuitive idea that points close to p are mapped by / into points 
closed to f{p). Prom equn.(13),(14) we have, 

X T 

\dt\sch= ^j^2 _r2)l/2 l^^l^'S'' M^l-gcfe = ~ (J^2 _ y2)l/2 I^^IkS (18) 
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and, 

X T 

\dr\sch = ^^2_y2)l/2 l^^l^'S' \Msch = ~ (X2_T2)l/2 l^^l^g (1^) 

where | | denotes the norm in the respective coordinate systems and we find that the coordinate transfor- 
mation, (t, r) — > (T, X) is not continuous on the horizon as the muhiphcative factors diverge on the horizon 
{X = ±T). Consequently the coordinate transformation {t,r) — >■ {T,X) is not a homeomorphism and the 
two coordinate systems do not topologically represent the same space-time manifolds [3,17]. Hence we show 
that that the Kruskal-Szckcrs coordinate system is not a proper extension of the Schwarzschild coordinate 
system and it is not a solution of the Einsteins equation for spherically symmetric black hole. We conclude 
this discussion with the following note: 

For any coordinate system we have. 

Consequently it is not possible to find a coordinate system with a regular g' without absorbing the 
singularities of {gsch.)piT at r = 2M into the transformation coefficients at r = 2M i.e, without breaking 
the diffeomorphic equivalence of the two coordinate systems. Thus, as also discussed in the preceding sections, 
the Kruskal-Szekeres coordinate system with a regular metric at the horizon can not be diffeomorphically 
equivalent to the Schwarzschild coordinate system and thus do not represent a static asymptotically flat 
solution of the Einsteins equations representing a black hole formed out of the gravitational collapse of an 
uncharged spherically symmetric asymptotically flat star, [see also AppendixrE] . 

It is not obvious how to describe the space-time evolution of the complete gravitational collapse of matter 
flelds as a whole in terms of time-like curves as, for a Schwarzschild observer, the time-like curves suffer a 
discontinuity across the horizon and become space-like inside the black hole event horizon. It is welknown 
that expressed in terms of the Schwarzschild coordinates the black hole event horizon has profound impact 
on the quantum description of matter flelds and black hole evaporation through Hawking radiation makes 
the space-time dynamic. Also Hamiltonian evolution of matter flelds break down on the flxed point sets of 
the time-like Killing vector field [14]. The canonically conjugate momentums are not well-defined on the 
horizon as will be evident from the lagrangians of the matter fields. 

V. DISCUSSION 

In a gravitational collapse once the collapsing body crosses the horizon it collapses to form the space-time 
singularity breaking the description of space-time in terms of continuous manifolds and the local symmetries. 
Wc can only characterize the presence of of the space-time singularity in a diffcomorphism invariant way, in 
terms of the curvature invariants along the space-time curves which cross the event horizon and necessarily 
terminate along the space-time singularity. The formation of black hole event horizon can be characterized 
through the formation of trapped surfaces. The gravitational collapse and the cosmological evolution are the 
only two processes in nature through which a three dimensional physical system collapses to zero dimension 
(forming the space-time singularity). Here through zero dimension we mean a point or a collection of points. 
We will illustrate this aspect in Appendix:B. 

Einsteins equations break down at the space-time singularity. This is something similar to electrodynamics. 
We can determine the electric field for a point charge using the Maxwells equations. But the field strength 
diverges and classical electrodynamics break down at the point charge (the corresponding quantum theory 
QED is not a resolution to this problem. It has its troubles associated with the point-like interaction terms. 
However experimental observations confirm that all the elementary particles arc of finite volume). The 
formation of space-time singularity is associated with finite volume to zero dimension transition for the 
corresponding collapsing body and the richest structure that we can attribute to zero dimension is that of 
an analogue of (compact) three dimensional generalization of the Cantor set [5,6] provided we generalize the 
description of the collapsing matter fleld through a proper quantum theory [a generalization of the Pauli 
exclusion principle]. 

There are two ways that one can reach zero dimension from flnite volume breaking the continuous topology 
of space-time manifold. One is through the point contraction mapping which requires an inflnite number of 
iterations which, together with the discussions in Appendix:B, is in accordance with the fact that time is 
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a continuous parameter. The other one is through the formation of an analogue of the Cantor set (or any 
other discrete manifold with different cardinality) in which case the underlying physical process to achieve 
zero dimension may be discontinuous. A discrete manifold may not always form a normed vector space, e.g, 
the set of points (n + x) on the real line, where n is an integer and a; is a fractional number, can not form a 
normed vector space as the difference between two points do not belong to the set. Also it is not physically 
obvious to talk of causal structure, defined through propagation of signals, in a discrete manifold unless the 
manifold is space- like and frozen in time (which is defined through physical processes). As discussed earlier 
according to the General Theory of Relativity charges associated with space-time transformation symmetries 
are global properties of a continuous space-time manifold as a whole whereas a manifold without continuous 
topology can only have space-time independent charge. 

To describe cosmological evolution and black hole evolution we will have to generalize and geometrize 
conventional quantum mechanics in a suitable way. In these respects the principal aspects to be critically 
studied, as will be discussed later in this article, towards unifying the general theory of relativity and 
conventional quantum mechanics is the concept of diffeomorphism invariance associated with the general 
theory of relativity and unitarity associated with the conventional quantum mechanics. 

The facts that the continuous topology of space-time break down at the space-time singularity [Ap- 
pendix:B] (indicating no well-defined observables associated with spatial transformations for the cosmolog- 
ically evolving or collapsing matter fields in the near zero dimension region) and that nature chooses a 
particular cosmology lead us to conclude that diffeomorphism invariance (which for large scale structure of 
space-time is equivalent to invariance under coordinate transformations) is not of so important (as it is for 
solar system microscopic physics) concern for the corresponding physical laws. Rather the fact Schwarzschild 
coordinates and the Kruscal- Szekers coordinates are not diffeomorphically equivalent indicates that an ap- 
propriate choice of a suitable coordinate system is most important. However we can express the generalized 
quantum theory in covariant form. This will help us to compare the generalized quantum theory with so- 
lar system quantum physics where the physical laws are invariant under inertial transformations and are 
formulated in a covariant way under the corresponding coordinate transformations. 

To generalize the conventional quantum mechanics we should take into account the following important 
aspects: 

(1) Special relativity made the concept of size for ordinary objects a relative one. The strong curvature 
effects near the space-time singularity will spoil the concept of dimension for the elementary particles forming 
the space-time singularity. 

(2) Quantum mechanics is the mechanics of quantum states which do not exist independent of their 
realizations at least in principle, i.e, through interactions with other quantum states. For solar system 
microscopic physics the fact that the elementary particles and bound states (e.g, atoms) formed by them 
are of finite volume has to be considered in the corresponding quantum state description as long as the 
continuous structure of the space-time manifold holds. 

(3) Every measurement process through state reduction is a non-unitary operation on the space of quantum 
states [7] . Near the space-time singularity the strong curvature will destroy the description of matter fields 
in terms of a unitary quantum theory. 

As far as the cosmological evolution is concerned (zero dimension to finite volume and finite volume to 
zero dimension transitions) no observer physics is the exact description of the evolution of the universe in 
the near zero volume region. A proper generalization of quantum mechanics may be non-unitary in the 
sense that the evolutions of the possible quantum states (if the space-time description of matter is given 
by a particular family of space-time curves representing possible particle trajectories [5]) representing the 
collapsing physical systems may be non-unitary. 

However, it is obvious that the collapsing physical system (which is a bound system through gravitational 
interaction) collapses to zero dimension violating the conventional quantum mechanics based on the uncer- 
tainity principle (e.g, the electrons in an atom obeying the quantum mechanical principles do not collapse 
on the positively charged nucleas) and follow the deterministic laws of general relativity. 

In the context of the above disciissions an important contradiction between the gcincral theory of relativity 
(describing the gravitational interaction) with conventional quantum mechanics (describing the microscopic 
interactions of the elementary particles) is the following: 

Positivity of the energy momentum stress-tensor together with the general theory of relativity leads to 
gravitational collapses [8] and space-time singularities [9] where a three dimensional physical system collapses 
to zero dimension (breaking the continuous space-time topology) whereas positivity of the energy-momentum 
tensor together with the canonical commutation relations lead to the Pauli exclusion principle (unless one 
introduces additional structures about the space-time singularity). 

We, living beings, are characterized by the fact that we can control some terrestial processes. But we 
can neither change the physical laws nor the cosmological evolution. Many descriptions we had made are 
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either approximations (unitary quantum mechanics, quantum statistical mechanics) or ideahzations (point 
particles for microscopic physics). The discussions in the preceding paragraph together with the facts 
that black holes contain no scalar hair [11], that there is no physical explanation of only recombination 
for the virtual particles (which "interpret" real effects as in the Casimir effect) to form loops in quantum 
field theory [Physically, in vacuum, even in Feynmann's summing over path scheme it is not obvious why 
particle pairs produced at one space-time point will only recombine at another space-time point. The 
otherwise should give abundance of particles and antiparticles. We will discuss this issue in some details in 
Appendix:C.] and that the universality of the minimum uncertainity relations are lost in the gravitational 
collapses and arc questionable in the solar system microscopic physics [12] lead to conclude that the proper 
avenue towards unifying these two theories and thereby explaining the cosmological evolution completely 
will be understanding the space-time singularity and extending the conventional quantum theory as the 
position-momentum Canonical commutation relations are in accordance with the corresponding minimum 
uncertainity relations [13]. 



VI. APPENDIX:A 



We can obtain the one-form {d(j))a from a zero-form (a scalar field) cj) in an explicit coordinate variables 
notation: 



(d'^)'^ = Ea|r[(dx)ar (21) 

here the range of the summation is the dimension of space-time and it does not represent an infinitesimal 
change in (j). When expressed in a particular coordinate basis [(dx)^]^ will be just dx'', a coordinate one- 
form, and the /i -th component of (d^)a is as an arbitrary tensor T, in its operator form, represented in 
a coordinate basis can be expressed as: 



, d ^, d 



.....(^)(^)-(dx'')(dx^)... (22) 
Here (gf^), (gf^-) are coordinate unit vectors and dx'*''' are coordinate one- forms. 

In terms of a coordinate basis the covariant d Alembaratian operator can be obtained from the invariant: 



(d«.(d«. = SE»-iJ^ (23, 

here the explicit summations are again over /x, v with the ranges same as above. When expressed explicitly 

in a coordinate basis the Lagrangian density of a massless scalar field is given by: gn^d'^4>d'^4) = g^^dfj_(i)d^(j). 
The infinitesimal change 5(f) of the scalar field (j) can be interpreted as the scalar product of the one form 
(d0)a and the infinitesimal vector line elements. 



VII. APPENDIX:B 



Quantum field theory is the quantum theory of fields. It gives the dynamics of fields, the quantum 
probability amplitudes of creation and annihilation of particles, in contrast to quantmu mechanics which gives 
the dynamics of the particles themselves obeying quantum principle. For the same boundary conditions these 
two descriptions match in the form of their kinematic solutions. Only for the free particle boundary condition 
the conventional interpretation of propagators in Q.F.T as giving the probability of particle propagation is 
in accordance with reality as the quantum probabilities are nowhere vanishing in both the theories. For 
microscopic particle physics experiments the free particle boundary condition is a good approximation in 
practice but ideally the field (f) (or the quantum mechanical wave function if)) is spatially confined within the 
experimental apparatus. In one dimension it is meaningless to say that a particle is propagating from one 
point to another if the probability of finding (or creation of) the particle is vanishing at some intermediate 
points. 

For loops in one dimension, the momentum-space calculations give the probability that a pair of particles 
with given four-momentums are created at one space-time point and a pair of particles are annihilated at 
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another space-time point with the same four-momentums. This feature is transparent if one consider ah 
possible space-time particle trajectories to form loops in one dimension which cannot be possible without 
the possible space-time particle-antiparticle (originated with given four-momentums and annihilated with 
the same four-momentums) trajectories crossing each other at least once in between any two given space-time 
points. Similar feature will be apparent if one interpret the loop as a particle encircling between any two 
given space-time points with the four-momentums at these two given space-time points remaining the same. 

Let us illustrate this feature following the standard literature. We consider situations where free particle 
approximation hold. In vacuum at a given space-time point x (in one dimension), the particle- antiparticle 
pair production probabilities with two-momcntums pi,p2 arc | exp(— zpi..xi)p and | cxp(— 2]32-2;2)P respec- 
tively apart from normalizing factors. Once produced the quantum mechanically allowed stationary state 
position-space wave functions that are available to the particle-antiparticle pairs are ipp{x) = exp(— ipi-ari) 
and i})Ap{x) = exp{—ip2.X2) respectively where x denotes space-time points. The quantum mechanical joint 
probability that the particles produced at Xi,X2 =0 with two-momentums pi and p2 can again coincide at 
a space-time point between x and x + dx is: 



P,(pi,P2) = iVi^ (24) 
here iVi is the relative pair creation probability at the space-time point a; = 0. 

Px{pi,P2) is independent of x,{pi,p2) and ^ as for free particle approximation L ^ oo although the total 
probability of coincidence is unity when integrated over all space-time points. Hence quantum mechanically, 
numerous amount of particle-antiparticle pairs should be observed in any microscopic experiment performed 
during finite time-interval if there would have been spontaneous pair creations in vacuum. 

In passing we note that a space-time formed out of loops in vacuum, closed time-like curves, as the source 
can not have an intrinsic time orientation in contrary to what is realized in nature. 

We next note that in non-relativistic quantum mechanics the total joint probability that two distinguishable 
particles with energy E and momentums (fc, —k; k = \n\ » 1) can coincide at some point x (here x is 
position only) is: 



P{aU) = J J ij sin^{kxi)S{xi — X2)sm'^ {kx2)dx2][dxif (25) 

here the integrals are performed over the interval [—L/2, L/2]. This expression turns out to be unphysical 
as the corresponding probability turns out to be unphysical [P{all) = 3/2]. Whereas classical mechanically 
the maximum value of P{all) can be nearly unity when the particles suffer impulsive elastic collisions to 
stick together and come to rest at some point a;o. 

Similar features as discussed in this section in the contexts of equ.(24) and equ.(25) will appear in three 
dimensions. 

In semiconductor physics, as charge carriers holes, are fictitious objects introduced for simplifications. In 
reality, quantum mechanically in a p-type semiconductor the motion of holes are out of the movements of 
the valence band or the acceptor level electrons. What is the proper explanation of the polarity of the Hall 
potential in a p-type semiconductor? The electrons experience forces {qvB) in the same direction as that of 
the holes and in terms of the motion of the electrons the polarity of the Hall voltage should be opposite. 



VIII. APPENDIX:C 



In this section we first consider the action for the gravitational field [20] : 



where Gab is Einstein's tensor. General Theory of Relativity intcrwinds incrtial mass (in general energy- 
momentum) of matter with space-time through the principle of equivalence and the dimension of the coupling 
constant k {k = 6.67 x 10~^cm^ — gm~^ — sec~'^) is completely determined in terms of only mass, space 
and time unlike, for example, in electrodynamics where one have another fundamental quantity (electric 
charge) to determine the dimension of the coupling constant. This feature is transparent if we compare the 
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Newtonian-limit of the general theory of relativity with the Coulomb's law. The consequence of this feature 

is the following: 

If we set c, ^ = 1 the dimension of k is length-squared ([/^]) and it is no longer possible to set fc = 1 as 
this will make the concept of space-time dimensions meaningless. Alternatively we could have set c,k=l 
(see footnote: page no. 269, [20]) then Planck's constant become dimensionful {[P]). 

However we can set Boltzmann constant (fc') = 1 by giving temperature the dimension of energy. In the 
reduced units c,w^,k' = l the gravitational action becomes: 



Sg = -Kj Gyf^ct^x (27) 

where K have dimension [Z^^]. 

We now conclude our discussion on Quantum Field Theory. In the rest of this section we will use the 
convention that < a\a> gives us probability density. 

We first consider non-relativistic quantum mechanics in one dimension. In position-space the normalized 
quantum mechanical wave function ij) gives us the probability amplitude. {ip*ip)dx gives us the probability 
of finding the particle within the infinitesimal length interval dx. For a free particle one adopts the delta 
function normalization scheme for the quantum mechanical wave function: 



f 

J — < 



V^^(a;)Vfe,(x)dx = 5(fci-fc2) (28) 



In this equation the left-hand side is dimensionless while the one-dimensional delta function has dimension 
of length [Z] as is obvious from it's definition: 



J — C 



f{k)5{k - ko)dk = f{ko) (29) 



for a regular fmiction f(k). 

It would be appropriate to replace the free-particle boundary condition by periodic boundary condition 
which is a reasonable approximation in situations where free-particle boundary conditions hold as for a large 
length interval the spacing between the adjacent values of the momentum allowed by the periodic boundary 
condition is negligible. 

In the reduced units (c, = 1) the action is dimensionless. The action for a complex scalar field is given 
by: 

1 = J Ld'^x (30) 
where the covariant lagrangian density for a massive field is given by, 

l|£^|^_m^ , (31) 
2^x^^x^' 2^ ^ ^ ^ 

Consequently the dimension of (f) {(f)*) should be inverse of length ([Z]~^). In the second quantization 

scheme < /3|0|q! > replaces the classical field [21] and the expression < q;|(/)*(/)|q; > gives the probability 
density of creation or annihilation of particles. For free-particle boundary conditions the Euclidean-space 
generating functional for a real scalar field is given by [22] : 



We[J] = NEe~^''^'^°''^^~i'^'' irili-a^8'^+m^+v"(M)iS(x-i-x-2)] (32) 

The terms in the logarithm giving quantum corrections are not dimensionless and the third term is not of 
the same dimension as of the first two terms. 

For a real scalar field confined within a finite volume box with periodic boundary condition and consistent 
with the second quantization scheme we have (equ.3.28,[10]): 

= yilEr . '^^"^'"- exp[g^(n,x^)] + -i^J] ^""^-^ exp[-g^(n,x/-)] (33) 
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in a relativistic theory a covariant normalization using four volume would be appropriate and the normal- 
izing factor should have dimension [Z^^]. Clearly the dimension of (j) ([^~^]) in equ.(33) is no longer 
and equ.(31) no longer gives us a lagrangian density. 

In other words the dimension of the scalar field </> as required from the action determining it's space-time 
evolution does not mach with the dimension required in the second quantization scheme in order that one 
can interpret < a|(^*(}()|a > as giving the probability density of creation or annihilation of particles. One can 
absorb the the normalizing factor into the the fock state by multiplying it by a factor with dimension 
as in the second quantization scheme < I3\4>\a > replaces the classical field. This will be in accordance with 
the probabilistic interpretation of the field ^ as we have, 



<a|<^*<^|a> = <a|(/.*|/3>^^</3|<^|a>. (34) 
where the sum is taken over all possible states. 

However this will violate the interpretation of the Fock states as quantum mechanically the normalization 
of probability density < a|a > in a relativistic theory requires that the each of the normalizing factors for 
the Fock states, whose number depends on the number of particles present in the Fock state, should have 
dimension 

We now consider fermions and electromagnetic fields. 

The covariant lagrangian density for each component of the free fermion fields (e.g, electrons-positrons) 
which is formed out of their causal space-time motion is given by: 

and dimension of ip is again {[l]^^). After linearization of the second order partial differential equation 
satisfied by tp we get the Dirac equation: 

{ij^d^ - m)^) = (36) 
Hereafter the following lagrangian density is used to study Q.E.D: 



L' = i^ii^d^ - e^A^ - m)V - ^F^^^JF^^i^) (37) 

with the dimensions of ip as determined above the dimension of the first part of L' is no longer that of 
a Lagrangian density and the action formed out of it is not dimcnsionlcss in the reduced units. Also the 
current density, eip^^ip (although the four divergence vanishes), do not involve any momentum operator and 
it is not obvious whether it is possible to have, in any approximation, the conventional interpretation of 
current density as charge-density times velocity from this expression. 

The Lagrangian density of a charged scalar field which is similar to the quadratic Lagrangian density for 
Q.E.D is given by, 

Lt = Lf + Lgm + Lint (38) 
This complete Lagrangian density for a charged scalar is gauge invariant only if we take Lint to be, 



Lint = -A^f + e^A^4,*cl) (39) 

the second term do not have a transparent interpretation unless we consider screening effects from classical 
electrodynamics similar to the corresponding discussions given in Appendix:H of this article. 
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IX. APPENDIX: D 

In this section we will illustrate the discussions in the context of equ.(20) in section IV. 
The metric of the two-sphere S'^{6, (j)) is given by 



ds^ = ^JTl^{dQ'^ + sin^ Qdc^^) (40) 
Here < ^ < tt and < (jf) < 27r. For the unit two-sphere we have, 

= d6|2 + sin2 6»#2 (4^^ 

with the ranges of 0, same as above. This coodinatc system have the following ill-behavedneses: 

(i) The coordinate </> suffers a discontinuity along some direction from 27r to 0. 

(ii) 4> is degenerate at the poles ^ = 0, tt. In spherical polar coordinate system (r, 0,<?!)) the point (r = 

c, = 0) where c is a finite constant is obtained through identifying any two arbitrary points on a circle 
characterized only through distinct values of </>. Similar construction is valid for the point (r — c, 9 — t:). This 
will be obvious if wc construct a two-sphere S'^{0, (f>) from a two-dimensional circular strip by identifying the 
inner-boundary and the outer-boundary to two distinct points [see the discussions below eqn.(50) regarding 
the reduction of eqn.(48) to eqn.(50) for points on the polar axis]. This construction can be generalized to 
higher dimension. 

(iii) The metric is singular at the poles [see the above discussions and App.B (a point can't be obtained 
from a one-dimensional line-element without breaking the corresponding continuous topology)]. 

When the S'^{9, (f) can be embedded in the three dimensional Euclidean space one can introduce Cartezian 
coordinate system (x, y, z) through the coordinate transformation: 

a: = sin 6'cos y = sinSsin^, z = cos^^. (42) 

here x^+y'^ + z'^ = 1. Although the metric is regular in the Cartezian coordinate system the transformation 
coefficients (^^r^) are singular at the poles and also at some isolated points on the x — y plane demonstrating 
the discussions below equ.(20). On the other hand to obtain the Spherical metric, singular at two isolated 
points, from the regular Cartezian metric one has to introduce coordinate transformation with singular 
transformation coefficients. 

We can also introduce two homcomorphic stereographic projections to coordinatizc 5^(0, (f) embedded in 
B? . The first one is from the North pole 6* = on the Equator plane to coordinatize the Southern hemisphere 
f < 61 < TT. We have, 

9 6 
X = cot-cos0, F = cot-sin^ (43) 

and this transformation is a homeomorphism sX 9 = -k. The second stereographic transformation is from 
the South pole ^ = tt on the Equator plane to coordinatize the Northern hemisphere < ^ < ^. We have, 

6 9 
U = tan -jcos (j), V = tan -sin ^ (44) 

and this transformation is a homeomorphism at 9 = 0. 

The transformation between the {X, Y) and {U, V) coordinate systems is a diffeomorphism at their common 
domain 9 = ^. The metric is also regular in these coordinate systems. 

However the transformation coefficients between {X, Y) and (9, <p) coordinates are singular at the South 
pole (^, ^ =0 at 9 = 7r),i.e, this transformation is not a diffeomorphism. Similarly the transformation 
(9, (f) — > {U, V) is not a diffeomorphism. It can also be shown that the transformations are not homeomor- 
phism at the poles. Thus S'^ is not homeomorphic to the infinite plane (i?^). Also a geodesic (an element 
of a particular family of reparameterization invariant curves) on can never come close to intersect itself 
whereas geodesies on S'^ are closed curves. 

We now consider the Robertson- Walker cosmological model. The space-time metric in terms of comoving 
isotropic observers is: 
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ds'^ = -dr'^ + a2(T)[#2 + sinV(«^6'^ + sin^ Odf/)'^) (45) 

here 0<^/'<7r, O<0<7r and < (f) < 2tt. The constant-time spatial three surfaces are compact 
(topologically S^) and there is no four-dimcnsional spatial geometry available to embed The spatial 

metric is singular along the closed line elements 9 = 0, 9 = tt including the two point-poles i/) = 0,7r. 

The discussions in section IV [c!qu.(20)] and in this section show that metric singularities cannot be removed 
by diffeomorphically equivalent coordinate transformations. Thus the black hole and the cosmological metric 
singularities are unavoidable aspects of nature. 



X. APPENDIX: E 



In this section we consider the electrostatic potential of a polarized dielectric system. The electrostatic 
potential of a polarized dielectric system is given by, 



here P is the polarization vector of the dielectric material and R is the vector joining the infinitesimal 
volume element dv carrying a dipole moment Pdv to the point of observation. It's magnitude is given by 
eqns.(46). We can reduce equ.(55) to a simpler form consisting two terms: one from a bound surface charge 
density ab and another from a bound volume charge density pb, 



y=J_/ P:da_^r Y^Pdv ^^^^ 

'^'^^0 J surf ace ^ '^'^^0 J volume ^ 

Here ab = P.h, n is the normal to the surface of the material and pb = — V.P. 

The total volume charge density in presence of a polarized dielectric medium is given by: 



p = pb + pf+ Psb (48) 

where we include cjb = P-fi in the free volume charge density as psb through the introduction of a proper 
delta function. For example in the case of a dielectric sphere we have, 

Psb{r,0,cj>) = ab{e,cj>)5{r-rs) (49) 

Here r, are radial distance (not vectors) and the delta function have dimension inverse of length. 
We then have: 

eoV.^ = p = Pb+p'j (50) 
where p' ^ = pf + psb and the divergence of the electric displacement vector D is given by, 



V.Z3 = p' f (51) 

Everywhere apart from the surface of the dielectric we have p' ^ = pf and the above equation (51) matches 
with the conventional expression for the divergence of D: 



V.D = pf. (52) 

The effect of Psb should be taken into the boundary conditions for D . This will also have consequences 
to obtain the energy density of a given electrostatic configuration in presence of dielectric mediums [23] as. 
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W.D - V.Do = Psb (53) 

We conclude this section with a few comments regarding the electrostatic field energy in presence of 
dielectrics. 

The electrostatic field energy in presence of dielectric mediums can approximately be considered to consist 
of three parts [24]: 

Wtot = \j D.EdT E.EdT = Wfree + Wspring + Wbound (54) 

here e = eo(l + Xe)- We briefly explain the three terms considering the reahstic case of a dielectric filled 
charged parallel-plate capacitor: 

i) Wfree is the energy to charge the capacitor to produce the configuration with a given electric field. We 
can regain this energy if we discharge the capacitor by connecting the two plates through a conductor. 

ii) Wspring is the energy required to increase the atomic/molecular dipolc moments or to polarize the 
atoms/molecules depending on, respectively, whether the atoms/molecules have permanent dipole moments 
or not. This energy will be regained as heat when we discharge the capacitor. 

iii) Wbound is the energy required to polarize the dielectric as a whole. The dipole-dipole interaction energy 
for two dipoles with dipole moments pi and p2 and separated by r is: 



U = - o [Pi-P"2 - 3(pl.f)(p"2-r)] (55) 

47reo r'^ 

U is minimiim when the dipoles are antiparallel and maximum when the dipoles are parallel. Consequently 
for any statistically infinitesimal volume of the dielectric (i.e, volume elements which are very small com- 
pared to the dimension of system but large enough to contain sufficient number of atoms/molecules so that 
microscopic fluctuations can be approximately averaged to zero) the orientation of the atomic/molecular 
dipoles will be as isotropic as possible. To polarize the dielectric we have to orient the atomic/molecular 
dipoles in near-parallel configuration in a given direction and supply energy to increase the electrostatic 
energy of the dielectric. This energy, Wbound, will be regained as heat if we discharge the capacitor. 



XI. SUPPLEMENT:! 



We will now study the behavior of the covariant Klien-Gordon equation in the near horizon limit. 

We will first consider the spectrum of the covariant Klien-Gordon equation in the (3 -|- l)-dimensional 
constant curvature black hole background which contains a one dimensional fixed point set of the time-like 
Killing vector field. This black hole space-time was obtained by Prof. M. Bannados, Prof. R. B. Mann and 
Prof. J. D. E. Creighton through the identification of points along the orbits of a discrete subgroup of the 
isometry group of the anti-de Sitter space-time. They used a static coordinate system where the constant- 
time foliations become degenerate along a particular direction apart from the black hole event horizon 
giving a one-dimensional fixed point set of the time-like Killing vector. The metric in the Schwarzschild like 
coordinates is given by. 



ds^ = [de^ - sin2 e{dt/lf] + ^ + r2#2 (56) 

where f{r) = (■^— p^)^. These coordinates arc valid outside the horizon (r > rh) for < 6* < tt and 
< < 2n. It is clear that the constant-time foliation becomes degenerate along the direction 9 — and 
9 = TT giving to a one-dimensional fixed point set of the time-like Killing vector field. 

The covariant wave equation of a minimally coupled massive scalar field is given by, 

^d^iy/^g^'-'d,^) - m'i, = 0. (57) 
The solution of the angular equation is given by, 
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where x = cos0, /x = iE"/ and (i^ + /i) ^ an integer. Here F{—u, u + — n; ^-^) is the hypergeometric 
function. This solution is throughout the angular range < 9 < n. Consequently the energy spectrum 
is continuous with divergent density of states. 

Wc will now illustrate that the divergent density of states is a characteristic feature of the fixed point set 
of the time-like Killing vector field indicating the breakdown of the canonical formalism of the conventional 
quantum mechanics. 

We will illustrate this feature in the context of the Schwarzschild black hole which contains a two dimen- 
sional fixed point set (the event horizon) of the time-like Killing vector field. Since the Hawking radiation 
through which the non-unitary black hole evaporation takes place originates mostly from the near horizon 
region we will consider the behavior of the spectrum of the covariant K-G equation in the near horizon 
region. Since the space-time foliation is static we will consider the stationary states. We will consider the 
radial solution of the covariant K-G equation. For the Schwarzschild black hole the constant-time foliations 
become degenerate at the black hole event horizon and it is not obvious to impose any consistent boundary 
condition on the horizon although the radial solution can be obtained through the WKB approximation and 
the spectrum is obtained using the semi-classical quantization condition (Brick wall model). However the 
conventional interpretation of the boundary condition is not obvious for extreme charged black holes. 

To obtain the energy spectrum one impose the following boundary conditions on the wave function: 

= (59) 
for r = 2M + h. The internal energy is given by [29] : 

^=15^- 

Similar behavior is also obtained when one considers the behaviour of matter fields in the Taub-NUT space- 
time which contains a zero dimensional fixed point (in the Euclidean sector) of the time like Killing vector 
field. In this case the angular solution (similar to spin-spherical harmonics) satisfies only the regularity 
condition that the angular part of probability current density integrated over S"^ is finite. This angular 
solution is similar to the spin-spherical harmonics. 

The divergence of the density of states as /i ^ as discussed above is a characteristic aspect of both the 
black hole event horizon and the cosmological event horizon. 

We now make some comments regarding relativistic quantum mechanics similar to App:D. For relativis- 
tically covariant normalization of the quantum mechanical wave function (or each component of the wave 
function for spinors) of we have, 
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i)*i}d''x = 1 (61) 



This indicates that in the reduced units (c, ^ = 1) the dimension of V' is inverse length squared [Z ^]. 
While the lagrangian leading to the Klien-Gordon equation is given by: 

The action determining the space-time evolution of '>\> is dimensionless in the reduced units. This gives 
the dimension of -0 to be in contradiction to that ([Z^^]) obtained from the normalization condition. 

In passing we note that to have a consistent time orientation for any space-time manifold, which through 
the principle of equivalence is form of existence of matter fields, particles should follow a particular family 
of reparameterization invariant curves in the space-time manifold [5]. 
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XII. SUPPLEMENTII: A FEW COMMENTS ON DOUBLE SLIT INTERFERENCE 

In this section we will give a causal description of double slit interference experiments. We will assume that 
wc arc working in the Cartczian coordinates in flat space-time where in the Lorentz gauge the electromagnetic 
potentials and the fields are conjugate harmonics and the fields obey the wave equations. 

This analysis of interference patterns is based on the fact that the velocity of the causal part of electro- 
magnetic waves is finite [Supp.III] and the corresponding wave-fronts travel with finite velocity. For a purely 
action-at-a -distance theory the two superposing fields are at the same state at any field point as they reach 
simultaneously from the two sources. Thus there would be no redistribution of intensity and we wont get 
interference fringe pattern. When the wave-fronts are traveling with finite velocity the wave fronts that 
simultaneously reach at a screen point from the slits should have originated from the two slits at different 
instants of time. We will find that the causal description of interference experiments together with the 
universality of interference fringe patterns naturally lead us to conclude that the velocity of light in vacuum 
is independent of the velocity of the source and the velocity of light is same in all directions. The causal 
interpretation is also useful to determine the velocity of light from the interference experiments. We will de- 
scribe the corresponding situation with wave-groups. In particular we will give the mathematical conditions 
for maxima and minima for interference experiments with finite wave-groups. We will also discuss some 
conflicts between the physical description of the interference process for the wave-groups and the Fourrier 
transform description of the wave-groups. These discussions are significant in the context of the foundational 
aspects of non-relativistic and relativistic quantum mechanics. In particular we find that a super-luminal 
velocity like the phase-velocity of the relativistic wave-function can determine a physical phenomenon like 
the corresponding interference fringe pattern which in this situation determine the probability density of 
finding the particle on the screen and hence the probable space-time motions of the corresponding particle 
passing through the slits. 

We will follow the diagram [13F] of section 13.3 of [26] for the following discussions. 

In Figure. 13F [26] the wave fronts that simiiltaneoiisly reach at the screen point P from the slits Si and 
5*2 should have originated from the two slits at different instants of time. The wavefront from 5*2 should 
originate at time earlier than the corresponding interfering wave front that has originated at ^i. The 
phase difference is: S ~ LoAt = = ^^2^ where lu and A are the angular frequency and the wave-length 

of the light respectively, S2A is the path difference between two interfering waves. 

If we consider a field point 'P' at any instant of time 't' and with positive values of X and Y coordinates 
Ei{x, y, t) from the upper slit is interfering with E2 from the lower slit that is similar to Ei that had passed 
through 'P' at time At earlier. Thus the interfering waves from the slit Si and 5*2 is given by: 

E\ = Asm {cot — kx) and E2 = Asm [uj{t — At) — kx] = Asm [uit — k{x + 6x)] 

These expressions are consistent with the fact that at a; = 0, i.e, at the slits we are considering wavefronts 
El at time t' and wavefronts E2 at time t' — At'. 

Hereafter we can continue the conventional procedure discussed in [37] to obtain the fringe pattern. 
Similarly, to reach a particular screen point P simultaneously the wavefronts from different parts of an 
aperture should have originated from wavefronts at the aperture at different instants of time differing in 
phases. This gives a causal description of the diffraction patterns. These discussions lead us to the following 
conclusions: 

(i) We should note that, with the frequency of the light being constant, we have the following condition 
on At: LjAt = tt for the first interference minimum and similar will be the situation with other minima and 
maxima. Thus At(s) remain to be the same for all experimental set-ups. The universality of the interference 
fringe patterns (the universality of the path differences) indicates that the velocity of light in vacuum is 
independent of the velocity of the source and the velocity of light is same in all directions. A similar analysis 
is given in [37] in the context of Michelson-Morley experiment. 

(ii) The interference fringe pattern can be used to experimentally determine the velocity of light with a 
source with known angular frequency (w) from the expression for the phase difference 6 = oj^^. 

We now consider the corresponding situation with wave-groups. When the coherence-length is much larger 
than the separation between the slits and the screen it is easy to extend the arguments for monochromatic- 
light to show that we should be using the phase-velocity to determine the interference pattern. 

Now let us consider the situation when the separation between slits and the screen is large compared to 
the coherence-length or is of the order of the coherence-length. We assume that the shape of the wave- 
groups originating from the two slits remain the same at any instant of time although they are 
traversing in different directions. To illustrate, the wave-group can originate from atomic transitions and are 
traversing in vacuum. Let the wave-group from the slit Si, Ei, has arrived at the field point 'P' at a time t. 
If we want to have the the first minimum at 'P' the wave-group from the second slit, E2, should be arriving 
at 'P' at a time At later such that Ei at 'P' is ahead of E2 in phase by an amount equals to tt. In other 
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words the shape of Ei and E2 at 'P' at the instant t + At should differ by half the wave-length with Ei ahead 
of £'2- The required time-interval At is — Thus we can proceed similar to the earlier cases as we have 

for the phase- difference: u)At = ^Ax = kAx and we have the first minimum when Ax = ^. Similarly the 
second minimum occurs when Ei ahead of E2 by one and half wave-length from the beginning of Ei . Here 
we have assumed that the lengths of the complete waves within the wave-group remain the same throughout 
the wave-group. If this is not the case the angular frequency will be varying along the wave-group and the 
mathematical calculations will have to be generalized accordingly for different orders of maxima and minima. 
The physical conditions for maxima and minima remain to be the same. For a weakly dispersive medium the 
above physical arguments regarding the path differences remain valid but in this case both the the angular 
frequency and the phase velocity is varying along the complete waves with different wave-lengths within 
the wave-group. The value of the wave-length A in the above discussions should be taken at approximately 
at the time when they both reach the point 'P', i.e, t + At and when the wave-group is weakly dispersing 
the wave-length may be defined through averaging over the group and the length of the wave-group should 
be defined through time averaging. Similar considerations are required when we consider the higher order 
minima. The above physical arguments for interference maxima and minima remain to be valid even when 
the wave-group is strongly dispersive. In terms of the group velocity the above physical arguments can be 
put in the following mathematical form which remains to be valid even for strongly dispersing mediums: 

For positive values of Y -coordinates the n-th order interference is maxima determined by the following 
condition on the group velocity vectors of the two interfering wave-groups, 

[iv,i~v,2)t].c^) = Em 

In the above equation t = when the center of the incident wave-group passes through the slits and i is 
running from one to n. Depending on the relative value of the length of the wave-group and the separation 
of the slits there can be a region where the interference pattern is governed by the phase-velocity. The 
wave-lengths arc in general different and may depend on time for dispersive wavc-gro\ips. Here we only 
assumed that the dispersion does not depend on the direction of propagation. Similar discussions are valid 
for negative values of Y -coordinates with the suffices interchanged. 

The above discussions are also in accordance with the significance of the coherent soiuccs. We demonstrate 
this with the simplifying assumption that wave-group contains a few numbers of complete wave-lengths and 
we can define an average wave-length and the phase-velocities are not differing much from one another 
throughout the wave-group so that we can assign a single value to the angular frequency and hence to the 
time-period. This is the situation with the wave-groups discussed in the above paragraph. Let the coherence- 
length of the wave-group be L. The approximate number of average wave-lengths in the wave-groups at the 
point 'P' is n = Here T is the average time-period of the wave-group at the point 'P' (we will also 

have to consider a time averaging over the duration of the wave-group at 'P' to define L and A if the wave is 
weakly dispersing) . Following the discussions in the preceding paragraph the interval of time through which 
the front of E2 should differ from the front of Ei in arriving at 'P' so that it just touches the end of Ei at 
'P' is r = ^ = :;7^7i:;7- Here A is the average wave-length of the wave-group at the point 'P' (we will again 

P 9 P \ ^ 

have to consider a time averaging over the duration of the wave-group at 'P' if the wave is weakly dispersing 
as is considered to define the average time-period T). For a wave-group which encompasses a few number 
of complete waves with the wave-lengths are not differing largely from each other we can take Vp = ^ and 
we have r = — which is the coherence-time as is expected. 

One should note that the the physical arguments presented in this section determine the interference fringe 
pattern for a wave-group even with very few numbers of complete wavelengths although the corresponding 
Fourrier transform can include an wide range of monochromatic waves differing in wave-lengths and may 
also in respective phase velocities. If we consider the Fourrier transform the interference fringe pattern is 
determined by the principle of superposition and the phase velocities of the monochromatic waves. We give 
a comparative study of these two descriptions for a wave-group with few numbers (of the order of ten) of 
same complete wavelengths propagating in a non-dispersive medium. We also assume that the amplitude 
remains the same or varies little over the complete wave-group. To have the first interference minima (there 
can be negligible residual intensity associated with the uncanceled half wave-lengths one at the beginning 
of Ex and one at the end of E2) as determined from the elementary physical considerations discussed in 
this article the wavelengths of the Fourier transform components should be mostly odd integral (including 
one) dividends of a basic wave-length which is same as the wave-length of the complete waves within the 
wave-groups. The spectrum gets even more restricted if we consider higher order minima. 

On the other hand the Fourier transform of a sin function centered at the origin is given by the following 
expression: 

a{k) = {^)sin{^)L - (^)sin(^)L 
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Here the wave group is given by sin(Q!x) extending from — {L/2) to {L/2). Thus the mathematical analysis 
is not in accordance with the elementary physical description althoiigh the physical description is what 
happens in reality. For only very long wave-groups for which mathematically the phase- velocity can be given 
by a sharply peaked value these two descriptions agree with each other. Thus light wave groups of finite 
extent in vacuum are fundamental in nature similar to soliton solutions. We will illustrate the above aspects 
later. 

For electromagnetic wave-group or other wave-groups like the ripple wave-groups the energy density is 
related with the value of the fields and the phase-velocity gives us the velocity with which energy is trans- 
mitted within the wave-group. When we consider matter-waves the probability density does not depend on 
the wave function but on the amplitude and phase-velocity within a wave-group does not have the corre- 
sponding interpretation as is mentioned above for electromagnetic wave- groups or ripple wave- groups. For 
matter wave-groups the particle flux is related to the particle velocity. However if the electron is described 
by a long wave-group, i.e, the electron is having a sharply defined momentum or the separation between the 
slits and the screen is small compared to the length of the wave-group the discussions at the beginning of 
this section regarding nearly monochromatic waves remain to be valid to describe the corresponding inter- 
ference fringe pattern. In these situations matter-wave phase velocity determines a physical phenomenon 
like the corresponding interference fringe pattern as is illustrated at the beginning of this section. These 
discussions are significant for matter waves as the phase-velocities are different for the non-relativistic and 
the relativistic wave- functions. The phase velocity of the wave- function for a relativistic electron is ^ [45]. 
In the case of Davission-Germer experiment with relativistic electrons, if the separation between the sources 
and the screen is small compared to thc^ length of the wave-groups the interference pattern is determined by 
the phase velocity which is greater than light and not by the group velocity as is conventionally referred [45] 
which is same as the velocity of the particle. The discussions made in point(i) in the context of the velocity 
of light are important as the angular frequencies are different for the non-relativistic and relativistic electron 
wave-functions (and consequently the corresponding time-intervals At for the first maximum will also be 
different). As discussed in the point(ii) we have, say for the first maximum, 5 = 2t:). Here 5*2^ is 

the experimentally determined path-difference, ui is the relativistic expression of the angular frequency of the 
electron wave-function. Thus we have a way to determine the super-luminal phase-velocity experimentally 
which determines a physical phenomenon like the interference pattern determining the probable space-time 
motion of the particle. The corresponding wave-lengths are also different and the experimental interfer- 
ence pattern should be in accordance with the relativistic expression of the matter-wave wavc-lcngth for 
the electrons. To illustrate these aspect, the condition for first interference maximum lead to the following 
relation between the relativistic phase velocity (vpr) and the non-relativistic phase velocity (vpnr) (with 
the assumption that the path-difference remaining the same in both cases): = ^=j. 

Whereas De-Brogle hypothesis lead to the following relationship between the two phase velocities: -^^^ = 
Here /? = f where v is the velocity of the electron. Thus the interference fringe pattern should be different 
with the non-relativistic and relativistic wave-functions. An experimental verification of the interference 
fringe pattern with a super-luminal phase-velocity provides a fundamental basis for the relativistic mat- 
ter wave-function. The otherwise will indicate that second-quantization is a more consistent description of 
the relativistic dynamics of microscopic particles. Similar situations arise in the context of the collapse of 
the quantum mechanical wave- functions and the EPR paradox [46]. In the interference experiments the 
wave-function is completely delocalized between the different maxima prior to collapsing to the point where 
the electron is striking the screen. Thus the collapse of the wave-function is a non-local phenomenon. In 
EPR paradox the the quantum mechanical wave-function predicts that the spin or polarization state of one 
of a pair of particles can change non-locally due to the corresponding change in the spin or polarization 
state of the other particle. The probability density of finding a particle at some point is determined by 
the corresponding matter wave-function signifying the physical reality of the wave-function. In this article 
we find that under certain circumstances the probability density and hence probable space-time motions 
of a relativistic particle is determined by the super-luminal phase-velocity of the corresponding relativistic 
wave- function. 

To conclude, in this article we considered a causal description of interference which is different from 
Huygene's description. We found that the causal description of interference experiments indicates that 
the velocity of light in vacuum is independent of the velocity of the source and the velocity of light is 
same in all directions. The causal interpretation is also useful to determine the velocity of light from the 
interference experiments. We outlined the corresponding situation with wave-groups an have given the 
precise conditions for maxima and minima for interference experiments with finite wave-groups. We have 
also also discussed some conflicts between the physical description of the interference process for the wave- 
groups and the Fourrier transform description of the wave-groups. These discussions are important for the 
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foundational aspects of non-relativistic and rclativistic quantum mechanics. We also find that a super-luminal 
velocity like the phase- velocity of the rclativistic wave- function can determine a physical phenomenon like the 
corresponding interference fringe pattern determining the probability density of finding the corresponding 
particle in the region behind the slits and hence the probable space-time motions of the particle passing 
through the slits. 

The above arguments indicates that the proper expression for fringe-shift in the Michelson-Morley ex- 
periments [26] is 2ojd^, where u is the frequency of light. We can not follow the procedure followed in 
[26] to obtain the path-difference as the velocity of light, as is assumed, is different in different directions 
and as discussed above it is not proper to calculate the phase-difference between the interfering waves by 
multiplying the path-difference by ^ directly. The fact that we can have the interference fringe system 
by allowing one photon to emerge from the source at each instant indicate that we have to introduce the 
corresponding electromagnetic wave description as fundamental and this description do not depend on the 
width of the slits apart from diffraction effects. 

For double slit interference experiments with electrons, the wave fronts are the position-space wave func- 
tions -0 and are of the form iVe*'^'^ '"""*). The interfering states are the states of the electron at the slits at 
two different instants of time. Only for the central maximum the two states at the two slits are the same. 

il) is the complete microscopic description of the electrons in the electron -beam interference experiments. 
It is so as there is no underlying ensemble representation of the probabilistic interpretation of ip. If we have 
to assume that, unlike the electromagnetic interference experiments, ip is vanishing at the slits for slit width 
less than that of the diameter of the electrons then the finite volume property should be a fundamental 
aspect in the microscopic description of the electrons. In other words ip gives a microscopic description of 
the electron upto a certain length. We can also illustrate this aspect with the following question: 

What is the; quantum mechanical description of a Radon atom in a rigid box when the distance of consec- 
utive nodes and antinodes of tp is equal to or less than the diameter of the atom? 

XIII. SUPPLEMENTIII: COMMENTS ON CLASSICAL ELECTRODYNAMICS 

In this Supplement we will make a few comments regarding the basic laws of Classical Electrodynamics. 

In section:A we will consider the multipole expansion of the Cartezian coordinate inverse distance elec- 
trostatic potential for a regular and extended charge distribution. We will show that a proper binomial 
expansion in terms of the spherical polar coordinates of the Cartezian coordinate inverse distance electro- 
static potential of a point charge not at the origin lead us to a multipole expansion in terms of the Legendre 
polynomials. The binomial expansion considered in this section is the most general and differs from the 
conventional expansion and gives well-defined result everywhere apart from two points. We will find that 
the potential expressed in terms of the Legendre polynomials is not well-defined at r = — r' and differs 
from the corresponding value obtained from the (-^) Coulomb potential. This leads us to conclude that the 
Cartezian coordinate (■^) Coulomb law for a point charge not at the origin and the expansion, equation 
(83), in terms of the Legendre polynomials are the same object only in a restricted sense. This will also be 
apparent from the discussions in Section:B. Nevertheless following the experimental observations both in the 
static situations as well as in the dynamic situations we have adopted the solutions expressed in terms of the 
Legendre polynomials to describe the potentials in the spherical polar coordinates. The discussions in this 
article lead us to conclude that the force law in the Cartezian coordinates, given by the inverse square law, 
and the force law in the spherical polar coordinates, given in terms of the Legendre polynomials, are different. 
The potential expressed in terms of the Legendre polynomials remains to be well-defined everywhere for a 
regular source. 

In section:B we will consider the Poisson's equation in the static situation. It is not possible to apply 
the integral version of the divergence theorem if the vector field is not well-behaved. We will demonstrate 
in this section that in the spherical polar coordinates an expression analogous to the integral version of the 
divergence theorem for a point charge in Electrostatics is valid when the point charge is at the origin. 

We will discuss the derivation of Maxwell's equation from the force laws. We find that in fiat space 
time it is in general not possible to derive the local Maxwell's equations in the familiar form from the 
inverse square force laws if we use arbitrary curvilinear coordinates. The equations for both the potentials 
and the fields contains non-local terms irrespective of the choice; of gauge. This is the situation with the 
spherical polar coordinates in flat space-time. In the spherical polar coordinates (r, 0, </>), the divergence of 
the electric field and the curl of the magnetic field contain non-local terms if the fields are expressed in terms 
of the corresponding inverse-square laws. Similar discussions are valid in the dynamical situations. In the 
dynamical situations the well-known Green's function, eqn.(120), does not satisfy both the inhomogeneous 
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and the homogeneous wave equations in terms of the spherical polar coordinates (r, 0, (/)). We will discuss this 
issue in Scction:IV. We will in particular consider the fields of an oscillating electric dipole. We will find that 
both the potentials and the fields obtained from the Green's function (120) do not obey the wave equations 
for points outside the sources and do not contain (i) radiation fields observed in the experiments. We will 
also find that the corresponding potentials violate the Lorentz gauge condition. This may also remain valid 
in curved space-time. 

In this article we give a formulation of electrodynamics in an arbitrary curvilinear coordinates for arbitrary 
sources with the potentials as fundamental variables on the basis of the following assumptions: 

Following the experimental observations, we adopt the principle of covariance to write the equations for 
the potentials in curvilinear coordinates and in curved space-time with reference to the physical laws that 
photons propagates along null geodesies in vacuum. 

The potentials satisfy the homogeneous wave equations in vacuum and the modal solutions of the wave 
equations creates or annihilates the photons in a second quantization scheme in the curvilinear coordinate 
system. 

The electromagnetic fields are related to the potentials by the conventional relations. However the elec- 
tromagnetic fields need not to satisfy the wave equations in vacuum. In the Cartezian coordinates both the 
potentials and the fields obey the wave equations in vacuum. This is because of the fact that in terms of the 
fiat space-time Cartezian coordinates the partial derivatives in the curl and in the wave equation commute. 
This may not be the situation in an arbitrary curvilinear coordinates where the corresponding partial deriva- 
tives may not commute. Mathematically the above discussions are related to the existence of conjugate pair 
of harmonics in an arbitrary curvilinear coordinate system in four dimensions. In the Cartezian coordinates 
where conjugate pair of harmonics exist the potentials satisfying the wave equations obey the Lorentz gauge 
and the fields also satisfy the wave eqiiations in vacuum. This may not be the sitiiation for the potentials 
satisfying the wave equations in an arbitrary curvilinear coordinate system where conjugate pair of harmon- 
ics may not exist and the corresponding electromagnetic fields may not satisfy the wave equations. This is 
the situation with the spherical polar coordinates in flat space-time. In the spherical polar coordinates only 
the divergence of B and Faradey's law remain to be valid for time dependent sources. These laws define 
the fields in terms of the potentials and should remain valid in a potential formulation of electrodynamics 
irrespective of the nature of the sources. The electromagnetic fields give the forces acting on a charged par- 
ticle. The homogeneous wave equations for the potentials reduce to the corresponding Laplace's equations 
in the static limit of static sources and c — >■ oo and the corresponding solutions expressed in terms of the 
Legendre polynomials describe the experimentally observed static potentials. We will later demonstrate that 
in terms of an arbitrary curvilinear coordinate system it is in general not possible to find a gauge so that 
the potentials satisfy both the wave equations and the Lorentz gauge condition. 

For static sources all the Maxwell's equations remain valid for points outside the source both in the 
Cartezian coordinates as well as in the spherical polar coordinates. In the static limit of time-independent 
sources and c — )■ oo the wave-equations for the potentials reduce to the corresponding homogeneous Laplace's 
equations. The current density should be steady, time independent and satisfy the corresponding regularity 
conditions discussed below equation (86). 

We now consider the expressions for the potentials within the source for both the static as well as the 
dynamic situations. In some situations the potentials within the source are given by the inhomogeneous 
equations. This is the situation in the flat space-time with the Cartezian coordinates and the flelds are given 
by the inverse square laws. On the other hand in some situations for points both inside and outside the 
sources the potentials are obtained from the homogeneous equations in terms of the multipole coefficients. 
The homogeneous equations do not contain the source charge-current densities. The relationship between the 
source and the fields are expressed by the multipole coefficients. The multipole coefficients are determined 
either from the assumed nature of the source or from the boundary conditions. The boundaries, the boundary 
conditions, permittivity and permeability will be different for points inside and outside the source. This is 
the situation with the spherical polar coordinates in flat space-time and the potentials are given by multipole 
expansions in terms of the Legendre polynomials. We will discuss the solutions in eqns. (83,84,127) which 
agree with the experimental observations and satisfy the homogeneous equations. Mathematically, the 
distributions for the potentials in the Cartezian coordinates are either the Green's functions of the Poison's 
equation in the static situations or the Greens functions of the four dimensional wave equations in the 
dynamical situations. On the other hand in the if we use the spherical polar coordinates, taking divergence 
of w.r.t r does not give us the same result if we take the divergence w.r.t R or r' (with a negative sign). 
This is clear from eqn.(73) and (78) and in terms of the spherical polar coordinates {r,9,(j)) the inverse- 
distance potential of a point charge not at the origin is not the Green's function of the Poisson's equation. 
This also follows from the discussions in Section:A where we have proved that and the corresponding 
expansion in terms of the Legendre polynomials, eqn.(83), do not agree with each other for all points. In 
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this article, in the spherical polar coordinates, the distributions relating the potentials to the source for both 
the static and the dynamic cases are constructed out of the modal solutions of the homogeneous Laplace's 
equation or the four dimensional homogeneous wave equations respectively. Only for the particular case of a 
static point charge at the origin the electric field satisfies Poisson's equation in (r, 9, 0) with the charge density 
given by a delta function at the origin and the field is given by the monopole term. For Maxwell's equations 
in curved space-time, the discussions in this article demonstrate that the diffeomorphism invariance of the 
inhomogeneous Maxwell's equations may become a non-trivial issue even if one have the inhomogeneous 
Maxwell's equations in a particular curved space-time. Similar as above, photons travel along mill geodesies 
in vacuum is a more fundamental principle in this respect. 

The discussions in this article are significant for the corresponding Lagrangian formulations and quantum 
field theory. In terms of the spherical polar coordinates the lagrangian does not contain the interaction 
term explicitly. In a second quantization scheme solutions of the form eqns. (127,128,129..) are similar 
to the interaction action in the Cartczian coordinates and gives the coupling between the source and the 
electromagnetic field. The stress tensor for the electromagnetic fields in the spherical polar coordinates 
is to be obtained using the Lagrangian formulation. We discuss this issue for time dependent sources in 
Section:IV. The stress tensor in this situation is formulated in terms of the potentials. 

We will conclude Section:B with a realistic description of the electrostatic behavior of conductors with 
reference to departure from the inverse square law. We will also improve the proof of the curl of B law and 
the no-work law for the electric field. We have also discussed the forms of the Maxwell's equations in the 
spherical polar coordinates, (r, 6, (p), with sources depending explicitly on the fields like sources obeying the 
Ohm's law in this section. 

In Supplement :X we will discuss the electromagnetic potentials in the spherical polar coordinates for time 

dependent sources. In this article we only consider accelerated charges. We will construct the solution 
for the potentials following the static situations. The potentials satisfy the four dimensional homogeneous 
wave equation in the spherical polar coordinates. For points outside the source the solutions are given by 
a set of new polynomials. We will give a few polynomials for field points outside the source. We will in 
particular consider the potentials and the fields due to an oscillating electric dipole. To the leading order 
the electromagnetic fields agree with the conventional expressions although the; power radiated differ from 
the conventional expressions [24] through a different dependence on the polar angle. For an oscillating point 
charge the conventional Green's function in the spherical polar coordinates, equation (120), docs not give 
(i) radiation fields in contrary to the experiments. Wc will also give the solutions of the homogeneous wave 
equation for the potentials for r < r' , where r is the field coordinate and r' is the source coordinate. The 
solutions reduce to the well-known multipole potentials expressed in terms of the Legedre polynomials in 
the static limit. The expression of the time-averaged power radiated from an oscillating electric dipole as is 
obtained in this article, eqn.(140), agrees with the observations. 



In this section we will consider the binomial expansion of the (■^) Coulomb potential in terms of the 
spherical polar coordinates (r, 0, (j)) in a way different from the conventional procedure. We will find that 
this binomial expansion of the Cartezian coordinate {^) potential of a point charge not at the origin in 
terms of the spherical polar coordinates gives well-defined result for all points apart from two points. 

We consider the electrostatic potential of an extended but localized charged distribution. We will consider 
points outside the source. Let r'{r' ,0' , (j)') be the position vector for an infinitesimal volume clement dv' 
within the source which makes an angle 9' with the positive Z polar axis and an azimuthal angle (j>' w.r.t the 
positive X axis. Let r(r, 9, (j)) be the position vector of the point of observation (P) making an angle 9 with 
the polar axis and an angle (j) with the positive X axis. The magnitude of the position vector R between dv' 
and P is then given by: 



A. Multipole expansion and the Legendre polynomials 



2rr'cos(7) 



(63) 



where, 



cos(7) — cos 9' cos 9 + sin 9' sin 9 cos{(t) — </>') 



(64) 



The electrostatic potential at P is given by, 



26 



^ ^J_f p{r',0',cl>'y^ sin 9' dr'dO' deb' 

In the conventional procedures [24] one expands ( 1 + z) ~ 3 in terms of z where z = — 2xcos 7 and x = 
If we consider points near the South pole (w.r.t the field point), the conventional binomial expansion is not 
valid in general for all r' . The absolute value of the ratio of two successive terms is "^^(^ z and the maximum 
value of z is 3. 

We now consider an alternate binomial expansion for the potential which is valid for all r' > r and lead 
to the Legendre polynomials. For any point outside the source we have the following expression for the 
potential: 



V{P) = -^l Pir')r"sinedr'd9'd^' ^^^^ 

We first expand [1 — 2j^pJ^^^ 003(7)]"^/^ in terms of 2 ^^2^^^/^^ cos(7). We obtain the following expression 
for the potential: 



T r 11 9rr' 3 9rr' 

We thereafter multiply the resultant series by (r^ +r'^) ^ and lastly we expand different powers of 
(r^ + r'^) in terms of ^. For points outside the source r > r'. This expansion is valid for all points r' < r 
apart from r = ztf*. For r' = r the potential is diverging and the corresponding binomial expansion may 
not remain well-defined. We will later consider the problem associated with r = —r'. 

The resulting expression gives us a series in terms of the Legendre polynomials: 



(68) 



Hereafter we can use the composition rule [23] for the Legendre polynomials to express the potential in 
terms of the Legendre polynomials in the field coordinates. 

We note that for f= —r' the above expansion is not well-defined. This is because the expansion of 
is not well-defined for a; = 1. To illustrate, let us consider the potential of a point charge, q, situated at 
z = a. The potential expanded in terms of Legendre polynomials for r > a is given by: 

nrl = ^^E^P^i'o^(^) (69) 

For z = —a the series in terms of the Legendre polynomials is indeterminate as is expected from the 
Leibnitz test [46]. It gives either zero or 

47reoa 

Whereas Coulomb's law give us the following expression for electrostatic potential a.t z = —a: 



Clearly the potential obtained at z = —a from expansion in terms of Legendre polynomials is not well- 
defined and is not in accordance with Coulomb's law expressed in terms of the Cartezian coordinates. This 
is related with the fact that the expansion of is not well-defined for 2; = 1 and equations (70,71) show 
that and the corresponding binomial expansion may not be the same object. (-^) is not the generating 
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function for the expression (83) for all points. We will later find that the inverse square law and, in general, the 
expansions in terms of the Legendre polynomials are solutions of different differential equations in different 
coordinate systems and it is again expected that they may not be the same object without restrictions. In 
terms of the spherical polar coordinates (r, 9, (j)) the inverse square law and the corresponding expansion in 
terms of the Legendre polynomials satisfy different differential equations. We will find in Section:B that in 
the spherical polar coordinates (r, 9,(j)) we can not interchange r with R or with r' (at the cost of a minus 
sign) while taking the divergence of (-^) and hence in the spherical polar coordinates {r,d,(f)) the inverse 
square law does not satisfy the Green's function equation for a point charge not at the origin. On the other 
hand with equation (83), for both r > r' and r' > r, we can interchange r with r' while taking the divergence 
of the corresponding electric field. Both the divergences arc vanishing and equation(83), the potential for a 
point charge not at the origin and expressed in terms of the Legendre polynomials, also does not satisfy the 
Green's function equation for a point charge not at the origin and is a solution of the homogeneous Laplace's 
equation. Only the (i) Coulomb potential for a point charge at the origin satisfy the corresponding Green's 
function equation in terms of the spherical polar coordinates (r, 9, (/>). We will discuss this issue in detail in 
Section:B. 

Following the experimental observations we continue to describe the potentials in the spherical polar 
coordinates in terms of the Legendre polynomials irrespective of the fact that in general these expressions 
satisfy the homogeneous Laplace's equation in terms of (r, 9, (f>). Also in the dynamic situation the potential 
for points outside the source is given by equ.(127), which in a second quantization scheme is consistent 
with the experimental fact that photons travel along the null geodesies in vacuum, and expression (127) 
reduces to the corresponding part of r > r' from eqn.(21) in the static limit of static source and c — > oo. 
For a point charge not at the origin the potential expressed in terms of the Legendre polynomials is not 
well-defined at f = r' and r = —r'. However for a non-singular volume charge density for which the 
potential -prrrPiicos -f)p{r')dv' and the corresponding electric field vanishes in the limit the source 

volume Sv ^ 0, the potential and the field expressed in terms of the Legendre polynomials remains finite 
everywhere and the field is also well-behaved at all points. We will discuss this issue in detail below equation 
(85). 

B. Discussions on Poisson's equation 

In this section we will consider the Poisson's equation for the electric field. We first note that we can not 
apply the integral version of the divergence theorem if the vector field is divergent at any point inside the 
region of interest. The expression Ex{Ax) = Ex{r = 0) + Ax[^^]r=o is not defined for a point charge at the 
origin. In one dimension we can use the relation: J^^ [df / dx)dx = f{x2) — f{x\) only if f{x) is well-behaved 
for all x between and including the limits. In the following section we will show that in the spherical polar 
coordinates an explicit volume integration of the divergence of the electric field of a point charge gives us 
the total charge when the charge is at the origin. Also for a point charge at the origin the integral of the 
divergence of the electric field is vanishing when the volume of integration does not include the origin. This 
together with the fact that the total fiux of the electric field of a point charge at the origin over an S"^ 
centered at the origin is given by ^ lead us to us a relation for the electric field of a point charge at the 
origin analogous to the integral version of the divergence theorem. 

We will shortly find that when the point charge is not at the origin the electric field in the spherical polar 
coordinates (r, 9, (j)) satisfying local differential equations is not given by the inverse square law. However 
in the following discussions we will ultimately consider a point charge at the origin and we use the inverse 
square law in {r,9,(j)) for the electric field of a point charge at the position vector a. We integrate the 
divergence of the corresponding electric field over a sphere centered at the origin. The radius vector of the 
point charge is a and the radius of the sphere is greater than |a|. The integral of the angular derivatives 
vanishes. For example the ^-integral of the polar angular derivative term is given by, 

A-Ke^Ie = £d9-^sm9 J E^^:^^{9.r')dv' (72) 

and this vanishes. Similarly the ^-integral of the azimuthal derivative term vanishes out of uniqueness. 
Here one should first integrate the source integral which contains a delta function at r' = a. 
The r-integral of the radial derivative term gives us the following result. 
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[r^ + — 2racos7]3/2 J ^ [^2 _|_ _ 2racos7]3/2 ^'^^^ 

In the limit a = this gives the result 1^ = j^. Hence we have JVp.Edv = ^ when the point charge is 
at the origin. 

When the point charge is at the origin 47reo JVp-Edv is vanishing if the integrating volume does not 
include the origin. If we consider any anmilar region centered at the point charge at the origin, the volume 
integral of the radial derivative term of Vf-E is given by the following expression: 



/ 



qdn[C-^f - C-f] = (74) 
f2 ri 



Whereas the ^-integral of the polar angular derivative term is, 

'd0§-^^^^9j P^0.f')dv' = O (75) 

we also have r' = 0. The azimuthal part vanishes out of uniqueness. It is easy to show that the volume 
integral of S/p.E is vanishing for any arbitrary volume not including the point charge at the origin. 
Thus we have the well-known Poisson's equation for a point charge at the origin, 

V.-.(^) = ATT5^{f) (76) 

The above discussions together with the fact that the total flux of the electric field over a closed surface is 
^ allows us to write the corresponding integral version of Gauss's divergence theorem for the electric field 
of a point charge at the origin. 

We can derive a corresponding expression when the charge is not at the origin following the above procedure 
through replacing r and r' by: 

{R, R') = (f - ?, - P) (77) 

respectively and considering the explicit integrals as above over a sphere centered at r'. The source integral 

is in terms of r" — r' . 

Thus we have the following relation, 

Vfl.(;§) = 47r<53(^) (78) 

where R = f— r' . 

We will later discuss that in general taking divergence with respect to the Cartezian coordinates do not 
give us the same result if we take the divergence with respect to the spherical polar coordinates. However 
the two coordinate systems give us the same result for the vector field (;&■). The divergence in both the 
coordinates are zero although there is a source at the origin. Thus we can use eqn(76) and eqn(78) in the 
Cartezian coordinates also. The experimental verifications of flat space-time Quantum Electrodynamics in 
the Cartezian coordinates also lead us to these relations. 

We should note that if we use the spherical polar coordinates, taking divergence w.r.t r does not give us 
the same result if we take the divergence w.r.t R or r' (with a negative sign). This is clear from eqn.(73) and 
(78). It is a unique feature of flat space-time metric expressed in terms of the Cartezian coordinates that we 
can interchange f and r' (at the cost of a negative sign) while taking the divergence of . The metric and 
hence R and the scale factors within the divergence operator are not symmetric w.r.t the field coordinates 
and the source coordinates if we use, say, the spherical polar coordinates. To illustrate the divergence of 
w.r.t f when f and r' are coUinear is geven by — ^r^^^rpp- The divergence of the same vector field w.r.t r' is 
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given by p^^^^pyr- It is obvious that the above two divergences are not antisymmetric. In fact in terms of 
the spherical polar coordinates (r, 9, 0), the Poisson's equation for the inverse square field takes the following 
form: 

fp.E = tr + t0 + t^ (79) 

Here, 

f p(r') d ,r^(r — r').r, , , 

'r=r-y9-ri^^w (80) 



and 

te 

and 



_ f Pi?) d {r'.0)siu0 



_ fpir') d .(r'.0) , 

*^-y ^a^^-R^^'^^ ^^^^ 

For a non-singular volume charge density the terms tr,t0, are in general logarithmically divergent when 
the field point is within the source and give non-local terms when the field point is outside the source. Also 
it is not possible to derive the other Maxwell's equation from the inverse square force laws and using the 
spherical polar coordinates. The equations for both the potentials and the fields contain strong non-local 
terms both in the static as well as in the dynamic situations. This is contradictory in a second quantization 
scheme where the modal solutions of the potential equations, as operators, are expected to create photons 
traveling in null geodesies in vacuum. On the other hand it is very well verified by experiments that in the 
static situations the expansions in terms of the Legendre polynomials describes the fields expressed in terms 
of the spherical polar coordinates for both points inside and outside the source. The multipole expansions in 
terms of the Legendre polynomials (more generally the associated Legendre polynomials) in general satisfy 
the homogeneous Laplace's equation in the spherical polar coordinates {r,9,(j)). Only for a point charge 
at the origin the potential satisfy the Poisson's equation. As we have discussed earlier the field of a point 
charge not at the origin is not given by the inverse square law in {r,9,(j)). The field can be found from 
the assumption that the potential satisfies Laplace's equation, from the fact that the potential is symmetric 
under the interchange of the field coordinates and the source coordinates (it can be easily seen from the 
following expression for the potential that the electric field is antisymmetric under the interchange of the 
source point and the field point when both of them are on the positive Z-axis and this relation is valid in 
general for the following expression as long as it remains well-defined) and from the boundary conditions at 
infinity. We first consider the simple situation of the azimuthally symmetric field of a point charge on the 
polar axis. The potential satisfying the above conditions are given by: 



In the more general situations where the point charge is not on the polar axis cos 9 is replaced by cos 7. 
Only for the particular case of a point charge at the origin the potential is given only by the monopole term 
and satisfies the Poisson's equation. We can obtain the field for an extended source by generalizing the above 
expression. For a regular source the field inside the source is given by the following expression: 

P I P 

4neE{r) = -Vr ^[^pP;(cos7)] - vW ^[-I^pP;(cos7)] (84) 

JrKr' ^ J r^r' 



When the charge density is non-singular the field is well-defined everywhere inside the source. The diver- 
gence and curl of E is zero. Thus we find that in terms of the spherical polar coordinates, for points both 
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inside and outside the source, the static potential and the field is given in terms of the Legendre polynomials 
which satisfy the homogeneous equation, the Laplace's equation. One should note that if physically the above 
solution gives us the correct field, then through Helmholtz's theorem, in terms of the spherical polar coor- 
dinates, the Laplace's equation will be the correct equation for the electric field. The relationship between 
the solution and the source is obtained from the multipole coefficients defined over the complete source. In 
day-to-day problems one usually specifies the source and the multipole coefficients are obtained from the 
assumed form of the source density. In a more general situations the multipole coefficients are to be obtained 
from the boundary conditions. The boundaries, boundary conditions, permittivity and permeability will be 
different for points inside and outside the source. At the origin f, 6, (p arc vanishing and degenerate. For a 
regular charge distribution one can find the field by studying the behavior of the field at points near the 
origin. Inside the source the velocity of light will be determined by the permittivity and permeability of 
the medium. We should note that in the spherical polar coordinates the divergence of the electric field is 
vanishing although there exist a monopole term in the electric field with the total charge dependent on the 
field coordinates within the source. The multipole coefficients are determined by integrating over the source 
coordinates and while taking the derivatives of the fields the multipole coefficients should be considered as 
constants. The divergence of the monopole field give us a delta function at the origin and hence is vanishing 
at the field point when the field point is not at the origin. If the charge density is regular the divergence 
is vanishing even if the field point is at the origin as p6v ^ in the limit Sv ^ 0. Only when we have 
a point charge at the origin the divergence of the electric field is non-vanishing which is not the situation 
with a regular charge distribution. However the monopole term in eqn.(84) gives finite flux for a spherical 
surface with radius r and centered at the origin. This is consistent with the integral version of Gauss's law if 
we chose the fields and the closed surfaces with some restrictions. For a sphere lying completely within the 
source and centered at the origin the field throughout the region should be given by the multipole coefficients 
evaluated at the boundary two-sphere. The integral law is also valid for an annular region centered at the 
origin with both the boundaries within the source or one of the boundary inside the source. The field within 
the annular region is given by the multipole coefficients evaluated at one of the boundary two-spheres. In 
eqn.(22) the multipole coefficients depend parametrically on the field coordinates and the application of 
the integral version of Gauss's divergence law is not a straight-forward issue. For points outside the source 
Gauss's integral theorem is valid for any anmilar region surrounding the source. 

Lastly, if we consider V in terms of the Cartezian coordinates = and this is in general not valid 
in the spherical polar coordinates while taking the divergence of a vector field. In the Cartezian coordinates 
Vp.E — R-E and hence Vp.E — -^S{r — r') for the inverse square field of a point charge not at the 
origin. This is not valid, as is evident from eqns (73, 78), if we take the divergence using the spherical polar 
coordinates and these aspects are not considered in standard literature, e.g, equations: (6.58, 6.59) [23]. In 
fact the the inverse square law is not the solution of the following equation expressed in terms of (r, 9, cj)): 



f,.E = (85) 
eo 

and the solution of the above equation, if exists, is not physically significant. As we have discussed earlier 
only for a point charge at the origin the above Green's function equation has a solution given by V{r) = j:^^- 
For a point charge not at the origin the potential is given in terms of the Legendre polynomials and satisfies 
the homogeneous Laplace's equation which is a local differential equation. As we have mentioned below 
eqn.(71) the potential is well-defined for all points apart from r = ±r'. The potential given by the Legendre 
polynomials is singular at r = r'. The potential at f = — r' is not well-defined. However both the above 
problems are not physically significant in a quantum mechanical description of charge densities or if we 
take into account the finite volume of the elementary particles so that the charge densities are finite and 
satisfy the regularity conditions discussed below eqn.(71). For non-singular sources the force law in the 
spherical polar coordinates (r, 6, c6) is given by eqn.(84) and for a non-singular source this expression is well- 
behaved for all points. Equation (83) is not the solution of eqn.(85) and satisfies the homogeneous Laplace's 
equation although the total flux of the electric field is non- vanishing for any sphere with r > r'. This is not 
contradictory to the integral version of the Gauss's law. We should note that the electric fleld is singular at 
the point charge and as we have mentioned earlier, in general, it is not possible to apply Gauss's integral 
law to singular flelds. We can apply Gauss's integral law to the corresponding electric fleld for any region 
leaving out the point charge. In the spherical polar coordinates (r, 9, (j)) we have an integral version of the 
Gauss's divergence theorem for a point charge only when the point charge is at the origin. An analogous 
statement similar to the integral law is also valid for the higher order pure multipole fields like the dipole 
field although the fields are singular at the origin. This is due to the fact that the pure multipole fields are 
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regular everywhere apart from the origin and the integral law is valid for any annular region surrounding 
the origin. The divergences of the pure multipole fields are vanishing everywhere. The fields are regular 
everywhere apart from the origin and since the radial dependence is of the form 1 jr^ ; 7t > 3 the angular 
dependence are such that the integral theorem remains to be valid for any annular region centered at the 
origin. 

We conclude our discussions on Maxwell's equations for steady sources with a few discussions on the 
derivation of the curl of B law. In the derivation of V x B law we have a 'boundary' term [24] : 

iB = \ ^p.[^]dv' (86) 

We can not use Gauss's divergence theorem directly as ^^^^ is divergent at .R = 0. We break the integral 
into two parts: 

Here Sv is an arbitrarily small sphere surrounding r. For a regular current distribution the first term 

vanishes in the limit i? — > 0. In the second term ^^^^ is well-behaved every where within V — 6v. We can use 
the divergence theorem to transform the volume integral into two boundary surface integrals. For the inner 

boundary J ^^^^.ds{ is vanishing for a regular current distribution as the radius of Sv become arbitrarily 
small. For the outer boundary J.dS2 is zero as the current should be tangential at the boundary of the 
source. However this derivation is valid when the derivatives of J are well-behaved within the source apart 
from the boundary. Similar discussions are valid in the spherical polar coordinates with (-^) is replaced by 
the expansion in terms of the Legendre polynomials. As we have mentioned earlier we can interchange the 
variables r and r' at the cost of a minus sign while taking the gradient of the potential (83) and together with 
the required regularities of the source similar to those mentioned above we have V^.A = in the spherical 
polar coordinates. We should note that the above discussions are valid only for steady currents with the 
required regularities. 

We can prove the divergence of B law from the fact that the magnetic field is the curl of the vector potential 
as is obvious from the Biot-Savart's law. However in the spherical polar coordinates the magnetic vector 
potential should be expressed in terms of the Legendre polynomials. It can be shown from the multipole 
expansion of the magnetic vector potential that there exists no magnetic monopole field [24] . 

We now consider the no- work law for the electric field. For regular sources the arguments will be similar as 
above for the magnetic field so that we can apply Stoke's law to prove the no- work law for the electric field. 
For a point charge we can consider an annular region surrounding the point charge with the inner boundary 
being an infinitcsimally small circle centered at the point charge. The work done for this inner boundary 
vanishes out of spherical symmetry and we have the no work law for the outer closed line. In the spherical 
polar coordinates we can consider an additional infinitcsimally small circle centered at the south-pole w.r.t 
point charge. The work done for this circle vanishes in the limit the radius of this circle vanishes and we 
again have the no-work law for the electric field. 

In the above discussions in this section on electrodynamics we have considered the potentials and the fields 
with static sources. We have found that with the inverse square laws the electric and the magnetic potentials 
do not obey the corresponding Poison's equations in terms of the spherical polar coordinates. In terms of 
the spherical polar coordinates the Laplacian of (-^) do not contain a delta function and give non-local 
terms. Following experiments we have assumed that the potentials for a point charge not at the origin is 
given in terms of the Legendre polynomials, equation (83). Equation (83) satisfy the homogeneous Laplace's 
equation and (-^) is not the generating function for the expression (83) for all points. Similar discussions are 
valid for the magnetic vector potentials for steady currents satisfying the regularity conditions mentioned 
below equation (86). In the static situation of static charge density and static steady currents we have the 
following equations for the potential and the fields in the spherical polar coordinates: 

V2F = 0; ^"^1=0 

the fields are given by: E = —VV and B = V x A. In the spherical polar coordinates all the Maxwell's 

equations remain valid only in the homogencoiis form. When time varying steady currents are present we 
will have to consider Faradey's law. The potential formulation will be consistent with the Faradey's law if we 
redefine the fields in terms of the potentials in the following way: E = —VV — dfA and B = V x A. However 
Maxwell's first and fourth equations, even in the homogeneous form V.E = and V x _B = fiQe^dtE may not 
remain valid in the spherical polar coordinates for time varying steady currents which give electromagnetic 
radiation. We discuss this situation in the following paragraph. 
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We now consider the most general situation of non-steady currents. To discuss the Maxwell equations 
for the potentials and the fields we will consider the radiation from an oscillating point charge in terms 
of the spherical polar coordinates in Supplement:X. We first discuss the conventional treatment in terms 
of a claimed Green's function of the wave equation in the spherical polar coordinates. We will discuss in 
Supplement:X that the Green's function, eqn (120), is not the solution of the corresponding inhomogeneous 
wave equation. This is similar to the discussions in the static situations below eqn (79). The resulting 
expression for the potential and the fields do not obey the wave equations. Also the Lorcntz gauge condition 
is violated. Thus the corresponding potentials and the fields are not the conjugate pair of harmonics in 
the spherical polar coordinates. We resolve this problem with the help of a second quantization scheme 
and the fact that photons travel in null geodesies in vacuum. We assume that there exist a scalar electric 
potential and a magnetic vector potential obeying the homogeneous wave equations in the spherical polar 
coordinates. In a second quantization scheme the modal solutions of the above wave equations creates 
photons traveling along the null geodesies in vacuum. The potential for a general source is obtained from 
the distribution formed in terms of the modal solutions of the wave equations for the potential. This is 
similar to the corresponding static situation where the potentials in the spherical polar coordinates are given 
by the distributions formed out of the normal modes of the homogeneous Laplace's equation expressed in 
terms of the Legendre polynomials. We will give a few modal solutions of the wave equation in the spherical 
polar coordinates. The electromagnetic fields are related to the potentials by the relations similar to the 
above paragraph: E = —VV — dtA and B = V x A. If the fields thus obtained would have satisfied the 
wave equations we woiild have a conjugate pair of harmonics in the spherical polar coordinates and the 
Lorentz gauge would also have been satisfied. In the spherical polar coordinates in four dimensions this 
appears not to be the situation. The resulting fields do not obey the wave equations although in terms of 
a more fundamental description, namely the corresponding quantum field theory, nothing acausal happens. 
The photons created by the modal potentials travel with the velocity of light. The well-known Maxwell's 
first and fourth equations for non-steady currents are not valid in the spherical polar coordinates. The most 
appealing aspect of the above formulation for time dependent sources is that the equations for the potentials 
reduce to the corresponding static equations in the static limit as we have discussed in Supplement:X. 

Thus in flat space-time we can not obtain the potentials and the fields in terms of the spherical polar 
coordinates by changing the coordinates in the expressions for the potentials and the fields expressed in 
terms of the Cartezian coordinates. We should solve the equations for the potentials expressed in terms 
of the spherical polar coordinates to obtain the potentials and the fields in terms of the spherical polar 
coordinates. For the time independent sources we do not have the corresponding inhomogeneous equations 
in terms of the spherical polar coordinates. We write the corresponding homogeneous equations in terms of 
the spherical polar coordinates. We thereafter solve the homogeneous equations for both points inside and 
outside the source. This remains to be valid for source densities not involving the fields themselves unlike 
sources following the Ohm's law. For points inside the source the permittivity and the permeability will be 
different from that of the vacuum. The relationship between the source and the solution is obtained from the 
globally defined multipole coefficients. This is not similar to the corresponding situation in the Cartezian 
coordinates where the fields are given by the solutions of the corresponding inhomogeneous equations with 
the locally defined charge-current densities as sources. For time dependent steady currents or for non- 
steady currents we may not have Maxwell's first and fourth equation and adopt the potentials satisfying the 
homogenous wave equations as fimdamcntal physical variables. The fields are obtained from the potentials 
through the following relations: E = —VV — dtA and B = V x A. 

The above discussions may also be significant for field theory in curved space-time where the scale factors 
within the divergence operator the and the partial derivatives of R w.r.t r and r' may not have any definite 
symmetry. In these situations we follow the above procedures to describe the electromagnetic fields in 
terms of the potentials. Even if the electromagnetic field in a particular curved space-time is given by the 
inhomogeneous equations the discussions in this article demonstrate that the diffcomorphism invariance of 
the inhomogeneous equations is a non-trivial issue and photons travel along null geodesies in vacuum is a 
more fundamental principle in this respect. In passing we note that the electric potential outside the source 
of a spherically symmetric asymptotically flat charged static space time is given by the equation: ^ = 

Here 'Q' is the charge as is measured from the field at infinity and grr = hr"^ . 

The above discussions are also significant for constructing the Lagrangian of the electromagnetic fields for 
both curvilinear coordinates and in curved space-time. In particular the Lagrangian of the electromagnetic 
field in the spherical polar coordinates in the static situation is given by only the free field homogeneous 
part. In the dynamical situation the Lagrangian density in the spherical polar coordinates is expressed in 
terms of the four potentials: A^^ = (^,^4) and is given by the expression: —-^Aa^^A'^'^. The stress tensor 
is expressed in terms of the four potentials and is to be obtained from the Lagrangian. This is significant 
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in the context of the corresponding quantiim field theory. Although the lagrangian does not contain the 
interaction term explicitly, in a second quantization scheme, the solutions in the dynamical situations [cqns 
(127..)] are similar to the interaction action in the Cartezian coordinates and gives the coupling between 
the source and the electromagnetic field. Only in the special situation when the point charge is at the origin 
the expression remains to be the solution of the corresponding Poisson's equation expressed in terms of 
the spherical polar coordinates and the Lagrangian will contain the source density. 

We conclude this section with the following discussions on the electrostatic behavior of conductors. We 
should note that with an ideal surface charge density the electric field within an arbitrarily-shaped conductor 
may not vanish. This follows from Coulomb's law. We can consider a simple situation of a cube with sides 
of length '2a' and is centered at the origin. Firstly, the externally given charge should be distributed equally 
on the six surfaces and the surface charge densities on the surfaces should be the same function w.r.t their 
respective centers. Secondly, the charge distribution on the surfaces should be a symmetric function under 
reflections w.r.t the centers of the respective surfaces. The surface configuration is symmetric under reflection 
w.r.t the coordinate planes and the field together with the surface charge densities on the top or the bottom 
surfaces can not change if we just rotate the cube about the Z-axis by an angle tt. Let us now consider a 
point on the Z-ax.is with coordinates (0, 0, e) where - — > 0. If the field has to vanish at 'P' we will have 
the following condition on the surface charge density cr{u, v): J J f{u, v)dudv = -^J J {u^ + v^)f{u, v)dudv, 
where f{u, v) = — 2lHiE) — ^ and {u, v) are the orthogonal Cartesian coordinates on the surfaces w.r.t their 

centers. One can obtain the {u^ + v^) factor in the second integral by considering rotations about the X- 
and Y- axes by angle ^ together with the symmetries of the charge distributions on the surfaces. Clearly it 
is not possible to find a non-trivial surface charge density satisfying the above condition. 

We now give a more realistic explanation of the electrostatic behavior of a charged conductor. We con- 
sider an arbitrarily-shaped thin metallic shell. If we do not represent the charges at the outside surface 
by an ideal surface charge density (which is the situation in reality) there will be higher order multipole 
fields which are directed inwards for electrons. To thc^ Iciading order the multipole potential is of the form 
r J p{r')cos {^)dr' dO! . The exact form of the multipole potential will depend on the nature and the symmetry 
of the source. The corresponding multipole force will push the free electrons within the core towards the 
inner surface. However this will create a small positively charged core in between. This core will impart an 
attractive force on the electrons and ultimately an equilibrium condition will be reached. If the multipole 
force is attractive then also a layer of positive charges will be created near the surface as the electrons in this 
layer are attracted towards the boundary and again an equilibrium will reached. The field within the region 
enclosed by the shell will be almost vanishing apart from microscopic fluctuations near the outer surface of 
the shell. This whole distribiition of the charges can take place within a few atom thick layer at the outer 
surface of the shell. These discussions indicate that any small departure from the inverse-square nature may 
not be observed in experiments similar to those discussed in the introduction of [23] . Thus the electric field 
can have a term of the form a{r)^ where a{r) « 1 and n > 2 with the force within a conductor vanishing. 
In this case it is not possible, in general, to have a Poisson's equation. Also these correction terms will 
introduce non-locality. For points outside the sources we will have an equation different from the Laplace's 
equation. However these terms may be significant to explain charge screening. 

We should also note that Legendre polynomials are either symmetric or antisymmetric around 6 = ^ 
and for a charge distribution on an arbitrary shaped conductor which can not be expressed as a sum of a 
symmetric and an antisymmetric part we can not expand the corresponding potential in terms of Legendre 
polynomials. 

We can construct a well behaved current density for the elementary charge carriers: J ~ p{[r — rc{i)\)v'c{t)- 
Here we have considered the charge density of the elementary charge carriers to be spherically symmetric 
and well-behaved, vanishing smoothly to zero around the boundary. The particle is a sphere of radius a and 
we can have a local charge conservation law leading to Maxwell's equations. 

Faradey's law, with the following convention (Section:?. 1.3 [24]), 

for motional emf or induced electric field the direction of the current or the electric field along a closed 
loop and the orientation of the enclosed surface giving the magnetic flux are related by the right-hand thumb 
rule, 

together with the above discussions and the current density equation (differential version of the electric 
charge conservation law) reproduces Maxwell's laws of Classical Electrodynamics. 

We conclude our discussions on the Maxwell's equations with a few comments on the forms of the Maxwell's 
equations in the spherical polar coordinates, {r,9,(j)), with sources depending explicitly on the fields. We 
consider a source obeying the Ohm's law, J = aE. In this case V.E has to be non-vanishing out of the 
charge conservation law and is given by: V.^ = - i|f. In this case the fourth equation in the spherical polar 
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coordinates, V x _B = leads to contradiction if we take the divergence of both sides. Since J depends 

expUcitly on the field E we follow the corresponding inhomogeneous equation and the fourth equation, in 
terms of the spherical polar coordinates, is given by: 

V X B = iiaE + ^e|| 

This leads to the following expression for the charge density if we take the divergence of both sides: 

p{t) = p{0)eM--J) 

and the first equation becomes, 

V J = ip(0)exp(-f 

in accordance with the corresponding equation in the Cartezian coordinates [24] . 

Wc now make some comments on ideal surface and line current densities. Let us first consider a surface 
with an wedge, i.e, two surfaces joined along a curve making an angle which can, in general, vary along the 
curve. If a surface current density originates at one surface the current can not propagate to the second 
surface as the only way that the current can have vanishing perpendicular components to both the surfaces 
is to propagate along the wedge (geometrically for two surfaces making an wedge the only way that the ideal 
surface current density can remain tangential to both the surfaces at the wedge is to flow along the wedge). 
The above arguments demonstrates that we can not have surface currents out of point charges flowing along 
an arbitrarily shaped surface as an arbitrary surface can be expressed as a collection of infinitesimal flat 
surfaces with non-parallel normals. Similar arguments as above demonstrates that we can not have currents 
in an arbitrary curve out of motion of freely flowing point charge carriers unless one applies forces everywhere 
tangential to the surfaces/curves. 

After that wc have reestablished Maxwell's equations wc should note that the finite volume of the elemen- 
tary charge carriers indicates that ideal line/surface charge densities cannot exist in nature unless the charge 
carriers can move with velocity c in one or two directions. Also to have ideal line/surface current densities 
the charge carriers should move with speed greater than c. Point charges can not exist in nature and the 
electromagnetic self-energies of the elementary particles are not infinite. 

Wc will now consider the following problem: 

Wc can put a point charge Q ~ —^q ai the center of a square with four point charges q at the vertices. 
The system will be in electrostatic equilibrium. Similar static equilibrium can be defined with three coUincar 
point charges g, g, — g/4 with —(7/4 midway. It can be easily shown that for any other moving observer 
(moving with velocity V w.r.t the rest frame of the five charges) this system will not be in equilibrium (the 

electric forces wont balance each other and there will also be magnetic forces) and the point charges will be in 

2 

relative motion with respect to each other. The magnetic forces being smaller by a factor can not balance 
the electric forces. Thus in the moving frame the conflguration is not in equilibrium and the charges will 
either fly away to infinity or collapse making the system non-existent in the moving frame. Similar situation 
will arise with three charges if one moves in a non-coUinear direction. However the above discussions are 
not significant in the microscopic regime because of quantum mechanical uncertainity principle and charge 
screening. 

Wc conclude this section with a few comments on magnetic monopolcs in the Cartezian coordinate system. 
It is straight forward to show that under a general electromagnetic duality transformation, eqn.(6.151, 6.152) 
[23], 

E = E'cos ^ + i?'sin ^ D = D'cos ^ + B'sin ^ 
H = -E'siuS, + H'cos£, D = -i5'sin^ -|- #cos$ 
and 

pe = PeCOS £, + p^sin ^ Je = J^cos ^ + J,'„sin ^ 

Pm = -PeSin^ + PmCOS^ Jm = "^Sm^ + 4n^0S$, 

the differential version of the magnetic charge conservation law doe's not remain time-reversal symmetric 

due to the pseudoscalar and pscudovcctor nature of magnetic monopole charge and magnetic current density 

vector respectively, i.e, if V.J^ = — is time-reversal symmetric then V .Jm = no-longer remains 

time-reversal symmetric although the differential version of the electric charge conservation law remains 
time-reversal symmetric under electromagnetic duality transformation. 



XIV. SUPPLEMENTIV: COMMENTS ON HYDRODYNAMICS 



In this article we will review the laws of fluid dynamics. Our discussions will be based on mainly that of 

chapter 40, 41 of The Feynman Lectures on Physics, Vol.2 [25]. 
The dynamics of dry water is governed by eqn.(40.6) [25]: 
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|^ + («^-VK=-^-V<A (87) 
or using a vector analysis identity to the second term of the above equation: 

^ + iVxv)xv = -^-Vcl> (88) 
at p 

where v is the velocity of an fluid element for which such laws can be applicable, p is the fluid pressure and 
(p is the potential per unit mass for any potential force present. We can derive some important laws from 
eqn.(86). The first one is the equation for vorticity (f2 = V x ■?) and is obtained by taking curl of eqn.(87): 

— + Vx(Oxt?) = (89) 
at 



The second one is Bernoulli's theorems (40.12) and (40.14) [26]: 



_^^,P , ^ , 1 2^ 



v.V{-+<j3+-v') = (90) 



i.e, 



— !-(/'+ -v^ = const {along streamlines) (91) 
P 2 

valid for steady flow and 

P 1 

— h + -v^ = const (everywhere) (92) 
P 2 

valid for steady and irrotational flow. 

However in all these equations the variation of the fluid density, p, is not considered while in deriving 
eqn.(40.17) [25] the variation of fluid density is not properly taken into account. The consideration for 
variation of the density of a nearly-incompressible fluid may become important through the facts that when 
unconstrained the shape of a fluid can be changed almost freely and sparsed away and through the facts 

that layers of fluids can be very easily spread or detached away although these properties vary from fluid to 
fluid. These features together with the local version of the conservation of mass law (assuming that there is 
no local source or sink in the region of interest): 



V.{pv) = (93) 

indicate that we should consider the possibility for variation of p properly as we will illustrate later that 

some ideal models can cause a finite variation of p and in reality the description of the motion should be 
changed. While the divergence of v may become important in cases like Couette flow where the centrifugal 
forces imposes a flnite and may even be large divergence of v. 

We can derive a proper version of eqn.(86) by applying Newtons second law to the fluid momentum per 
unit volume and we have: 



+ [v.V] ipv) = -Vp - V{p<P) (94) 

This equation is in general a non-linear coupled [through eqn.(91)] partial differential equation for v. 
Bernoulli's theorems for fluid dynamics can only be established when p is constant : 
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i.e, 



v.V\p+{cl>p) + \{pv^)]=Q (95) 



p+ {(pp) + -{pv^) = const {along streamlines) (96) 



valid for steady flow and 

1 



p + {(j)p) + 2^pv^) = const {everywhere) (97) 
valid for steady and irrotational flow. 

In general, when p is varying, only the flrst of the Bernoulli's theorems : 

v.W\p+{<t>p) + \{pv^)] = i) (98) 

remains to be valid provided {v-Vp) is vanishing or is approximately valid if | {v.Vp)v\ is negligible compared 
to the other terms in eqn.(93). To illustrate the significance of these comments, let us consider the ideal 
model to calculate the efilux-coeflicient, fig. 40-7 [25]. After that the contraction of the cross-section of the 
emerging jet has stopped we have, from the conservation of mass law, the following equation for pv at two 
vertical points: 

PlVl = P2V2 (99) 

In this case pressure is the atmospheric pressure and remains the same throughout the flow and thus even 
for the flow of a nearly-incompressible fluid p can vary as v changes with height. In reality the flow usually 
gets sparsed away after a distance which varies for different flows. 

The viscous flow of a fluid is governed by the following two laws which are obtained from eqn.(93) and 
eqn.(41.15),[25]: 

+ [u.V](piT) = -Vp - V(p(^) + r]W^v + {r, + r)')V{V.v) (100) 



V.{pv) = (101) 

supplemented by proper boundary conditions. To illustrate the significance of the boundary conditions we 
can consider the change of the shape of the surface of water in a bucket when the bucket is given a steady 
rotational motion about it's axis. The surface of the water become paraboloidal when the bucket is rotating. 
This shape can not be obtained without a vertical component of fluid velocity along the bucket surface for 
a flnite duration although the bucket surface only have an angular velocity. 

In the above equations rj is the "first coefficient of viscosity" or the "shear viscosity coefficient" and rj' is 
the "second coefficient of viscosity" . This equation is extremely significant in the sence that this equation, 
not eqn.(41.16) [26], is the equation which contains all the terms relevant to describe the dynamics of viscous 
fluids, both nearly-incompressible and compressible. For compressible fluids p will also depend on pressure, 
p{r). We can modify this equation only through varying the nature of the viscous force. 

The equation for vorticity is given by: 

— + V X (O X tT) - n{V.v) + VX V(V.?T) = ^(V^fl) 
_ (VpWM _ [V(P^)i X (VP) ^ ^(V^^^) X (VP) ^ !L±!/ v(V..-) X (VP) (102) 
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We can obtain an equation similar to cqn.(41.17) [25] describing the motion of a viscous fluid past a 
cylinder provided we can neglect the terms involving Vp and it is given by: 

— + V X (j] X u) - n{v.v) + VX y{y.v) = ^(v^fi) (i03) 

Following the procedure in section 41-3, [26] we can rescale the variables to obtain an equation which has 
Reynolds number (i?) as the only free parameter : 

^ + V' X (tJ X u) - cI;(V'.u) +UX V'(V'.u) = ;^(V'^cD) (104) 

where the prime describe the scaled variables, w is the scaled velocity and R is given by the usual expression, 

R = ^VD. 

p 

To conclude in this section we have derived the exact equation describing fluid dynamics. We considered 
the motion of both non-viscous and viscous fluids. We proved that in both the cases there are terms which 
arc neglected in the conventional theory but may become significant in some ideal model and in reality 
the description of motion is changed. Some of these terms even change the dynamical laws of viscous fluid 
motions by violating the conventional theory established in term of the Reynold number and these terms 
are signiflcant for the dynamics of compressible fluids like air. 



XV. SUPPLEMENTV: A FEW QUESTIONS 



What happcins to the entropy increase principle as the Universe evolve to form the big-crunch singularity? 
What happens to the uncertainity relations along the process of gravitational collapses? What is the quantum 
mechanical description of a Radon atom in a rigid box when the distance of consecutive nodes and antinodes 
of ijj is equal to or less than the diameter of the atom? What is the position-space wave function of two finite 
volume massive bosons if we take contact interaction [we can characterize contact interaction in a relativistic 
covariant way if we take Ax, At = 0, i.e, do not consider the coordinates of the center of the particles 
to characterize the interaction kinemetically] into account? How a photon produce electron-positron pair 
with finite volume concentrate rest masses? What are the charges and masses of the electron-positron pairs 
forming loops in the vacuum? How two particles with three- momentums fci, ^2 (fci ^ ^2) produced to form 
a loop at a space-time point always arrive at another spacial point simultaneoTisly? What is the microscopic 
explanation in terms of particle exchanges of the force in the Casimir effect? (If a loop begine at one plate 
and terminate at the other then the pairs, being tangential, can not impart any energy-momentum to the 
plates. This is also relevant for the discussions at the beginning of this article and discussions regarding 
bound states in [38]) What is the mechanism of the collapse of the momentum-space wave function of a 
particle knocking out an electron from an atom? What is meant by j^* >= ci{t)\'i>(j2?,s > +C2(i) |5't/i234 > ? 
Is the radiation out of a single atomic transition a spherical wave? Quantum mechanically the region between 
the rigid walls (which is equiprobable in classical mechanics) is non-homogeneous for a particle in a rigid 
box ! A single photon can not reproduce Maxwell's equations apart from moving with velocity c. How can 
a process involving only a few photons be described starting from the Maxwell's equations? The large scale 
structure of the Universe is homogeneous. 



XVI. SUPPLEMENTVI: A FEW COMMENTS ON VIEWING 



Let viewing be a relative process which is achieved by comparing the images of the surroundings formed 
at the brain with the images formed at the brains reduced by a factor N. We can scale time- intervals 
without measuring any spatial length. Time-intervals are primarily and fundamentally determined by our 
brains through the axial rotation of the Earth and do not involve any change at our brain similar to the 
spatial extensions of objects while seeing. Consequently the velocities and accelerations of the images of 
particles and signals at our brains are reduced by a factor N. Thus the images (and their dynamics) of the 
surroundings at our brains are not an exact replica of the surroundings. Relating the surroundings with the 
corresponding images formed at our brain is a part of mind/inteligence. 

Is the images formed in our brain are always reduced in size? 
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One can adders to this question through elementary observations on babies. In this context we should 
also note that forming larger images is essential to a have developed distant vision. 

We note that with a sharp resolution (lO^^m if we consider that we can see a 0.1mm thick hair at Im 
away and the focal length of the eye-lens to 1mm) the processing of image at the retina can involve only one 
cell. 



XVII. SUPPLEMENTVII: A FEW COMMENTS ON UNCERTAINITY PRINCIPLE AND 
RADIATIONS INSIDE ISOTHERMAL ENCLOSURES 



In this section we first discuss the position-momentum uncertainity relation in the context of optical 

microscope observations. For microscope observations if we analyse the scattering of a photon by an electron 
through Compton scattering then to the leading order the uncertainity in the x-component of the electron 
momentum for photon wavelengths greater than twice the Compton wave-length of the electron is 



2/ICOS0 2h^ 1 ^ 
Ap^ = -^cos(^ (105 

Here is the angle of scattering of the photon in the forward direction and the photon is assumed to be 
incident along the positive X-direction (7r/2 — <;;!) is the angle of incidence of the scattered photon w.r.t the 
microscope lence) . The maximum uncertainity in of the electron, in the context of Compton scattering, 
that the electron is observed happens for </> = 0. 

Consequently the approximate minimum position-momentum uncertainity relation becomes: 



AxAp^ = h (106) 

mocA 



The expression (for </> = 0) valid for all A is: 



^ ^ 1.22/1 1.22/1 2Ac , , 

AxAp, = 107 

n 2n A -I- 2Ac 

Here Ac is the Compton wave-length of the electron. 

The r.h.s varies with A and for A >> Ac the right hand side of the minimum uncertainity product is 
. The above equation gives the minimum uncertainity product in the sense that we have considered the 
maximum resolution for the incident photon with wave-length A and this leads, in principle, to a value ^^^^ 
for the r.h.s when A << Ac- 

In general the right-hand side also depends on cos (p. 

We now discuss a few aspects of radiations inside isothermal enclosures. According to Wien's law [30] the 
radiation inside a hollow iron isothermal enclosure at it's red-heat temperature is peaked about a wave-length 
which is about five times larger than that of visible light and the radiation in the visible spectrum may be 
negligible. 

Are the radiations inside a slate-isothermal enclosure and an iron-isothermal enclosure the same at the 
red-heat temperature of iron. 

In the derivation of Raley- Jeans law for radiation inside a cubic isothermal enclosure we have the following 
relation between the allowed set of modes (with perfectly reflecting boundary conditions) for two cubic 
enclosure of sides L' and L respectively: 



Now if there exists a radiation within the second cube for a particular set of values of (ux, ny,nz) the above 
equation in general may not give integral set of solutions for {n'^,ny,n'^) for arbitrary values of L and L'. 
Thus the corresponding mode of radiation will be absent in the first cube, something which is unexpected 
from the radiation mechanism through atomic structure with the material remaining the same. 

The above discussions are also relevant for the equilibrium quantum statistical mechanical description, in 
the context of varying the volume of the enclosure, of particles in a rigid box. 
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We now consider an interesting paradox. We can consider two absorbers inside an isothermal enclosure. 
The absorbers can only absorb radiation with two different bands of wave-lengths and pass the rest. Since ex 
is different we can construct a Carnot-engine similar to those discussed in page-118 [30] violating the second 
law of thermodynamics. However in reality this do not lead to any contradiction. Let us consider the process 
in detail. Let the absorbers be parallel to two faces of the cube. Due to the presence of the absorbers there 
will be tubes of radiation with a given range of wave-length (the portions that the absorbers absorb) that 
wont be reaching corresponding area-elements of the parallel surfaces. But these area-elements will keep 
on radiating the corresponding radiations. This will remain valid at any stage. Thus the temperature of 
the enclosure will keep on decreasing and it is not possible to operate the Carnot-engine without changing 
anything else and thereby violate the second law of thermodynamics. These arguments are also applicable 
for the Carnot-engines considered in page-118 [30]. In this context he perfectly reflecting mirror models are 
misleading. 

We will later discuss a few aspects of the quantum statistical concept of entropy for an ideal gas kept 
inside an isolated enclosure taking the collision interactions into account. 

In passing we note that the expression TVe'^"^*"*^^) leading to negative energy do not satisfy the non- 
relativistic Schrodinger equation. We will discuss later the relevance of this aspect in the context of the 
convergence of the perturbation coefBcients with time for perturbing potentials sinusoidal in time. 

XVIII. SUPPLEMENT VIII: A FEW COMMENTS ON NEWTONIAN MECHANICS 

In this section we will make a few comments regarding Newton's laws of motion for conservative force 
fleld. 

We define a conservative force to be a force for which the infinitesimal work done by the force field is total 
derivative, i.e, 

F.dr = -dV (109) 

Here dr is an infinitesimal element on the trajectory of the particle in the given force field and is determined 
by the definition of force in terms of the velocity and its time-derivatives. 

The l.h.s of the above equation, F.df, will be a total derivative if F become either m-^{v) or 

The first choice gives us Newton's second law. 

Let us consider the terms of the second category. These terms do not take into account the vectorial nature 
of velocity completely. The force remains to be vanishing either if the velocity is zero or if the magnitude 
of the velocity remains constant. If the velocity only changes its direction then this definition of force leads 
us to conclude that there is no force acting on the particle. The simplest example is the rotational motion. 
Thus these definitions are not in accordance with the first law of motion. 

Another possible class of total integrals are of the form Z?^ [^ •''7] • Again the terms inside the parentheses 
are of the form ^ {y.v) and the force will again be vanishing if only the magnitude of the velocity remains 
constant. 

We wont consider terms involving time-derivative of the acceleration. The force will then be vanishing for 
constant accelerations. Also we wont consider r dependent terms as these terms will identically lead to zero 
force either at the origin or at infinity even if the acceleration is non- vanishing at the corresponding points. 

If one can show that the definitions of force as are considered in this section are the only possibilities so 
that F .dr is a total derivative and the first law remains valid, then we can say that the second law follows 
from the first law at least for conservative force fields. We should also note we can not define force to be a 
sum of two terms one of which satisfies the above criteria but the other does not. 

In passing we should note that the equation F = —VV{x, y, z) is valid only for the solutions of the above 
equation joining two points in a given interval of time;. Other trajectories give us different accelerations and 
thereby different forces acting on the particle. If we do not define force to be mass times the acceleration 
the definition of kinetic energy will be different. We will later consider the significance of these discussions 
in the context of the Lagrangian formulation of classical mechanics. 

We should note that piecewise discontinuous curves involves impulsive forces and one should take into 
account this aspect while discussing conservative forces. 
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XIX. SUPPLEMENT IX: A FEW COMMENTS ON THE ENTROPY OF A SCALAR FIELD IN 
SCHWARZSCHILD BLACK HOLE BACKGROUND AND THE BRICK WALL MODEL 

The thermodynamical nature of black holes was first established by Bekenstein [39]. He obtained an 
expression for the entropy of black holes. This expression was proportional to the area of the horizons of the 
black holes. Hawking obtained an exact expression for the entropy by considering the behavior of matter 
fields in the black hole backgrounds [40]. The entropy of a black hole, considered as a thermodynamical 
system, was found out to be Here A is the horizon surface area. The thermodynamical properties of 
matter fields in black hole backgrounds was discussed by 't Hooft [29]. The black hole is assumed to be in 
thermal equilibrium with surrounding matter. For a static black hole the rate of particles radiated by the 
black hole should be equal to the rate of absorption of matter by the black hole, 't Hooft [29] assumed that 
the matter field wave function is vanishing near the horizon as well as at infinity, i.e, for a large value of 
the radial coordinate. Since then a lot of works had been done on the entropy of matter fields in black hole 
backgrounds using the brick-wall model [41,42,43,44]. 

In this section we will discuss a few aspects on the entropy of a massless scalar field in the Schwarzschild 
black hole background using the brick- wall model of 't Hooft [29] . We will show that with a proper coimting 
of the density of states, the entropy of a massless scalar field is different than that obtained by 't Hooft. 

The radial part of the wave equation of a massive scalar field in the Schwarzschild black hole background 
is given by the following expression: 

(1 _ —)E^i,{r) + irdrr{r - 2M)dri>{r) - {'^^ + m'Uir) = (110) 
Here I is the angular momentum quantum number and E is the energy. 

We follow the procedure of t' Hooft to calculate the entropy of a scalar field using the brick wall model 
and the WKB quantization rule [29]. We consider massless fields but the discussions can also be extended 
for massive fields. The boundary conditions on the scalar field are the following: 



ij{2M + h),ij{L) (111) 

where h « 2M and L » 2M. Here h is the brick- wall cut-off parameter. We consider the solutions 
of the wave-equation obtained by the WKB approximation. The solutions are assumed to be of the form 
p{r)cxpiS{r). The solutions arc stationary w.r.t the radial variable throughout the spatial manifold. These 
solutions are also important for the analysis of the Hawking radiation. 

The expression for the free energy of a massless scalar field is given by the following expression: 



Here V{r) — 1 — The range of E integration is from zero to infinity although one should be considering 
the back reaction to the metric. The range of r integration is from 2M + h to L where L » 2M. 

We now consider the free-energy calculation using a proper density of states. Let us first consider the WKB 
quantization condition together with the brick wall boundary condition. We have the following condition: 



7rn= / drk{r,l,E) (113) 

Jri 

Here the symbols have their usual meaning. We assume that for a given set of values for E, I, ri and L 
there exist one solution of the above integral equation. The upper value of I is restricted by the condition 
that k{r) should remain real throughout the total region of integration. We have, 



MO = ^^i^^-^^a + l) (114) 

Here we are considering a massless scalar field. At r = 2M-\-h with h ^ 0, I can be varied freely. However 
for a state for which k{r, E, I) is real throughout the complete region from r = 2M + h to r = L the upper 
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limit of I is restricted by the maximum value of ^^^^ which occurs at r = 3M. Hence the upper limit of / is 
given by the following expression, 



1(1 + 1) = 27E^M^ (115) 

2 

The I integration in cqn.(llO) for the free energy introduces a factor yjyj in the radial integral. Wc can 
take the upper limit of / -integration to be l{l + 1) = s{27E'^M^), where s can be taken to be arbitrarily 
close to one but not one so that k{r) is real for all values of r. The product s^^^(27i^^M^) is less than 
for all values of r. We first do the / -integration. The free energy is given by the following expression: 

We can do a binomial expansion in s to find out the leading divergent term in the brick-wall cut-off 
parameter. The maximum value of the term 27M^^5^^ is one. We thereafter do the radial integration and 
the energy integration respectively. It can easily be seen the term associated with the linear order term in 
s is the most divergent term in the brick- wall cut-off parameter h. The entropy is given by 5 = /3^^. The 
entropy associated with the linear order term in 's' is given by, 



Where I = /o°° (c^pa:-i) ~ JE- "^^^ divergent part of the entropy associated with the brick- wall cut- 
off parameter is logarithmically divergent in the cut-off parameter. This is expected as the constant-time 
foliations intersect at the horizon. Here ^ is given by the usual expression /3 = SttM and h is the brick-wall 
cut-off parameter. Wc have assumed L >> 2M. Wc find that with this expression of /3 the divergent part 
of the entropy in terms of the brick-wall cut-off parameter is independent of the horizon surface area. The 
higher order terms in s are not divergent in h and vanish in the limit L — >■ oo. We can extend the above 
discussions to massive fields. For massive fields k{r) is given by the following expression: 



k{r) = jj^\lE^-mW{r)-^^l{l + l) (118) 



y(r) 



We first find out the maximum value of the term [yj^ + + 1)]- Thereafter we determine the 

maximum allowed value of l{l + 1) so that the square-root remains real for all values of r and I. The free 
energy is given by the following expression: 



-P^ = -l^j [eMm-lJ ^)[{E'-m'yir)}^-{E'-m'Vir)-s^NiN + l)U (119) 

Here A'' is the max;imum allowed value of I. We can now follow the steps similar for the massless fields 
with E"^ replaced by E"^ — m'^V{r). The leading divergence of the entropy in the brick wall cut-off parameter 
remains to be of the same form as that for the massless field. The energy is defined with respect to an 
asymptotic observer and the lower value of the limit of energy integral is m. 

To conclude, we have considered the entropy of a scalar field in a Schwarzschild black hole background. We 
have followed the procedure formulated by 't Hooft. The scalar field is vanishing near the horizon as well as 
far away from the horizon. The density of states is calculated using the WKB approximation. We found that 
the entropy associated with the states for which the WKB quantization rule is valid throughout the region 
of interest is different from that obtained by 't Hooft and as is discussed in the standard literature. The 
laws of black hole thermodynamics are differential laws. The entropy of the matter fields may be important 
to interpret the entropy renormalization constant [29] . The part of the entropy associated with the brick- 
wall cut-off parameter is independent of the black-hole mass and can be considered to give the entropy 
renormalization constant. However this term is divergent in the brick-wall cut-off parameter and leads to 
continuous spectrum. Also the universality of the part of the entropy associated with the brick- wall cut-off 



42 



parameter indicates the ill-behavedness of this part of the entropy is associated with the intersection of the 
costant-time fohations at the horizon as is mentioned earher in this section. In this context it is interesting 
to study the corresponding situations with other black holes backgrounds. The entropy of the matter fields 
is also important in the scattering matrix approach of 't Hooft where the entropy of the black hole may be 
related with the states of the matter fields outside the black hole [29] . This may be significant to explain the 
black hole entropy if we consider the black hole and the surrounding matter fields outside the horizon as a 
single system. However we have found in this article that the matter field entropy is not proportional to the 
horizon surface area, as was obtained in the earlier works, and the above interpretation is not valid for the 
entropy of the matter fields in the states considered by 't Hooft [29] and in this article. It is also important 
to note that the entropy of the matter fields in the Euclidian sector is proportional to the area of the horizon 
and differs from the Lorentzian sector expression obtained in this article . The expression of the matter field 
entropy for the states considered in this article is also important for the entanglement entropy approach to 
explain the black hole entropy [44] . The fact that the solutions of the matter field may not contain a global 
oscillatory term for all values of 'I' is also important to discuss the scattering and absorption of matter fields 
in black hole backgrounds. These discussions may also be significant in the context of quasinormal modes. 
We will consider the above aspects later. 

XX. SUPPLEMENT X: RADIATION FROM AN OSCILLATING ELECTRIC DIPOLE 

In this section we will consider the solution of the electromagnetic wave equations in the spherical polar 
coordinates for time varying sources. We will derive an expression for the the electric and magnetic potential 
of an oscillating point charge and an oscillating electric dipole in the spherical polar coordinates. 

In the conventional literature the potentials and the fields in the spherical polar coordinates are assumed 
to satisfy the non-homogeneous equations similar to the Cartezian coordinates. In the Lorentz gauge the po- 
tentials satisfy the non-homogeneous wave equations. The Green's function in the spherical polar coordinates 
described in reference [23], eqn 6.64, is given by the following expression: 

^/ ^ N 1 f expilcoT — kR] , 

G{R, = -'^ (120) 

Here t = t — t' . The uj -integration nms from — oo to +oo and k is always positive. However similar to 
the static situations, in terms of (r, 6, (p) the above equation does not satisfy either the homogeneous wave 
equation or the corresponding Green's function equation. It is easy to find that Vf?G — contains 
terms of the form J J{R) '['^"'^-'^1 (j,^ leading to no n- vanishing retarded terms in the wave equations. For 
regular sources the above Green's function leads to the well-known retarded potential: 

Vir,t) = j'-^^ (121) 

Where U = t— ^ \s the retarded time. We can also show from this expression that in terms of (r, 9, </>), the 
potential does not satisfy the either the homogeneous or the inhomogcneous wave equation thus violating 
the consistency of the above solution. Only in the Cartezian coordinate the retarded potential satisfies the 
inhomogeneous wave equation. 

We now construct a solution of the Maxwell's equations in the spherical polar coordinates. We follow 
the discussions made in the context of the static situations. The Maxwell's equations in the spherical polar 
coordinates defining the fields in terms of the potentials are given by the following homogeneous equations: 

y.B = (122) 

and 

fxE = -— (123) 

As we have mentioned in the introduction we here give a potential formulation of electrodynamics in 
the spherical polar coordinates with time dependent sources. We consider the potentials satisfying the 
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homogeneous wave equations as fundamental physical variables. This is consistent with a second quantization 
scheme and as wc will observe below give us the experimentally observed static and dynamic potentials and 
fields. If there exist conjugate pair of harmonics the fields will also satisfy the wave equations and the 
Lorentz gauge condition will be satisfied by the potentials. However in the spherical polar coordinates this 
appears not to be the situation. The potentials that we construct from equations (127,128,129..) do not in 
general satisfy the Lorentz gauge in the spherical polar coordinates and the fields obtained from equation 
(127,128,129..) may not obey the wave equations in the spherical polar coordinates. 
The homogeneous wave equation for the potentials is given by the following expression: 



1 

VV-^^=0 (124) 

The most simple solutions of the homogeneous wave equation in the spherical polar coordinates is given 
by the following expression: 



g{f,t) = NY,fiir,^)Pr{cosj)[expiojt] (125) 
Here Pj™(cos7) are the associated Legendre polynomials and fi{r,uj) satisfies the following equation. 



|:(r^f )-[U/ + l)-^r^]=0 (126) 

The construction of the distribution in the spherical polar coordinates is not a straight-forward issue. We 
have to consider the static distribution to construct the distribution in the dynamic case. Here we consider 
the Z = 0, 1, 2 terms for r > r' . fo{r,uj) is given by the well-known expression: 

Mr,.) = ao'-^^^tM (127) 
r 

/i(r, w) and /2(r, w) are given by the following expressions: 



/i(r,w) = ai[ik + -J 



(128) 



and 



. , , , ik 1 ,exp(— zA;r) /,„„s 

/2(r,u;) = a^i-- + - + ^] > (129) 

Here 01,02,... are constants. These constants are determined in terms of the sources and the resulting 
expressions give us the potentials in the dynamical situations. With the normal modes given by the above 
expressions we can determine the coefficients ai,a2, ... from the static limit of the above solutions. In the 
static limit the sources are time-independent and c ^ 00. In the static limit the above solutions reduce 
to expressions similar to the static multipole potentials if we choose the coefficients 01,02, ... properly as is 
demonstrated in the following discussions in this section. This is also the situation for points inside source, 
r <r' . Although the wave equation for the potential reduce to the Laplace's equation in the static limit the 
solutions of the wave equation may not always reduce to the solution of the Laplace's equation in the static 
limit. We will discuss this situation later and now consider sources for which the above solutions reduce to 
the corresponding static solutions in the static limit. For localized sources varying periodically with respect 
to time we define the coefficients 01,02, ... in terms of the static multipole coefficients determined from the 
amplitudes of the sinusoidally varying components. To illustrate, in the following we will be considering the 
radiation from an oscillating electric dipolc. In this case wc chose ao = q and ai = p where q is the charges 
of the oscillating point charges and p is the amplitude of the oscillating electric dipole moment. We also 
chose N = The angular frequencies of the radiation are determined from the source frequencies and 

the corresponding wave numbers are determined from the relation a; = cfc as is valid in the vacuum. 

The ^ electric field obtained from the above potentials in general does not satisfy the wave equation in 
the spherical polar coordinates. This is easy to obtain if we consider the angular dependence of the solutions 
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of the wave equations as are given above. Thus the potentials and the electromagnetic fields are in general 
not conjugate harmonics in the spherical polar coordinates. 

For points inside the source the solutions are given by different expressions and satisfy the homogeneous 
wave equations with velocity dependent on the permittivity and permeability of the medium. This is similar 
to the corresponding static situation in the spherical polar coordinates. We will later discuss the solutions 
of the homogeneous wave equation in the spherical polar coordinates. 

We will now consider the potentials of an oscillating point charge and an oscillating electric dipolc in the 
spherical polar coordinates for points outside the source. We have two point charges q, —q oscillating along 
the Z -axis with the coordinates given by the following expressions: 



5 S 

s+ = -cos{u)t')k; S- = —-cos{u)t')k (130) 

The dipole moment amplitude is p = qs. The charge densities are of the form p = q5^{r' — | cos(wt')A:). 
In the spherical polar coordinates the potentials of the oscillating point charges that are obtained from 
the Green's function (120) are given by expressions similar as below: 



47reo (l-^2)i|^_f?(i^)| 

Here (t^) is the retarded source time: = t——. However the dipole potential obtained from the Green's 
function (120) in the spherical polar coordinates does not give (i) radiation fields. On the other hand we 
will find that the distributions constructed in this article give us (i) radiation fields for an oscillating electric 
dipole. 

In the Cartczian coordinates the retarded potential due to the oscillating positive charge, obtained from 
the Lienard-Wiechert potential, is given by the following expression: 



V{r, t) = -^ ^-^ (132) 

47reo(l-/3.ii)|r-s+(ifl)| 

where tn = t — \J—I±ii^^ ^ = £. The potential for the negative charge is given by a similar expression. 
The complete expression in the above equation does not satisfy both the homogeneous and inhomogcneous 
wave equation in the spherical polar coordinates if we replace the Cartezian coordinates by the the spherical 
polar coordinates. This is similar to the corresponding electrostatic situation. However in the limit of small 
source velocities, in the dipole approximation A >> s and for {x,y,z) » {x',y',z') the total potentials to 
the linear order in l3 and (^) is similar to the corresponding expressions that are obtained from distributions 
(127,128,129..) if we replace the Cartezian coordinates by the spherical polar coordinates. The total radiation 
potential in terms the spherical polar coordinates and linear in the dipole moment is given in (equation 9.24 
[24] . It may be so that the Lienard-Wiechert potential for an oscillating point charge does not obey the wave 
equations. We will later discuss the derivation of the Green's function of the wave equation in the Cartezian 
coordinates given in [23] . It is more appropriate to discuss the radiation problems using the quantum theory 
of fields. 

In the spherical polar coordinates the relationship between the field coordinates and the source coordinates 
is given by equations (127,128,129..) and also by the higher order solutions. The potentials obey the wave 
equation by definition. 

We will first find the potentials obtained from the monopole term /o(r). The corresponding potential for 
an oscillating point charge is given by the following expression: 



Vo = Re\— exp(w< - kr)] (133) 

47reo r 

This term gives us the static monopole potential in the static limit. In the static limit we recover the 
Laplace's's equation valid for a point charge not at the origin. With an oscillating point charge the right 
hand side of the wave equation of the monopole potential vanishes for points outside the source r > s. For 
points r < s the distribution will be different. This is similar to the corresponding static potential in the 
spherical polar coordinates of a point charge not at the origin. The static monopole potential is vanishing 
for an oscillating electric dipole. We should note that if the point charge is oscillating with a sine function 
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then the corresponding monopole dynamical potential is vanishing in the static limit and the dynamical 
solutions do not reduce to the static solution. For a point charge oscillating with the cosine function in eqn 
(130) replaced by sine function this expected also from the fact that in the static limit the potential of a 
point charge at the origin satisfy the inhomogeneous Poisson's equation. However in this case also we define 
the coefficients ai,a2,... in terms of multipole moments determined from the amplitudes of the sinusoidal 
oscillations. 

We now consider (i) radiation potentials. These terms are obtained from the (i) terms present in the 

normal modes fi{r, k) with / > 1. The (i) radiation potential linear in the dipole moment is obtained from 
/i and is given by the following expression: 



y = i?e{^^exp[£l;i-fcr](cos6»)} (134) 

With the assumed form of the time varying sources the electric potential is given by the following expres- 
sion: 



V = _^(cos0)^5^^ii^l-^ (135) 
47reo c r 

Here k — kf and Cj = ck. The above expression is similar to the corresponding expression obtained in [24]. 
The second term in /i gives us the static dipole electric potential in the static limit of w — >■ and c — >■ oo. 

The (i) terms present in the higher order fc)s also give (i) radiation fields. However in the perfect 
dipole radiation limit, A >> s, the amplitudes of these terms are of higher order smallness and we neglect 
these contributions. 

The current density is given by J = pv. Each component of the magnetic vector potential satisfies the 
homogeneous wave equation in the spherical polar coordinates. We can obtain the magnetic vector potential 
following the above procedure with the coefficients ai,a2, ... determined in terms of the magnetostatic mul- 
tipole coefficients. To the leading order in the dipole moment p the radiation vector potential is obtained 
from fo{r) term and is given by the following expression: 



- ^0 . sin ((Di - fcr) . 

A = puj ^ -z (136) 

47r r 

The higher order terms in the distributions (127,128,129..) do not contribute radiation terms linear in p 
in the oscillating electric dipole magnetic potential. The above expression is similar to the corresponding 
expression obtained in [24]. 

We now consider the derivation of the electric and magnetic fields from the potentials. We take the 
derivatives of the potentials using the Cartezian coordinate unit vectors and transform the unit vectors in 
the resulting fields to the spherical polar unit vectors. The resulting electric and the magnetic fields, to the 

2 

order linear in (^^), are given by the following expressions: 



2 

E = -^^^E^(sme)sm(LJt - kr)e (137) 
47rr 



and, 



B = ^^-^—{sme)s\Q.{ujt-kr)^ (138) 

Here w is the angular frequency of the harmonic oscillator. The electromagnetic fields give the forces acting 
on charged particles and currents. However the electromagnetic fields in the spherical polar coordinates 
obtained from equations (127..) do not obey the wave equations and are not conjugate harmonics of the 
corresponding potentials. We will later discuss the polarization properties of the electromagnetic fields. 

We now consider the power emitted by the oscillating electric dipole. The Poynting's vector is to be 
obtained from the free field Lagrangian expressed expressed in terms of the four potentials = (^,A). 
The Lagrangian density is given by the expression: —-^Aa^jiAajj. The time averaged radial power obtained 
from the potentials (135) and (136) is given by the following expression [47]: 
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Thus the time-averaged radial power, the expression analogous to the radial Poynting vector, coincides 
with the conventional expression [24]: 



However the time-dependent energy-momentum tensor may be different from the conventional expressions 
[24]. This is expected from the fact that in general, in an arbitrary curvilinear coordinate system, the 
potentials may not satisfy both the wave equations and the Lorentz gauge condition. It will bo convenient to 
consider the potentials from the individual oscillating charges to verify that the potentials give a conserved 
energy-momentum four vector. We will later discuss this issue. Equation (140) is obtained by replacing the 
Cartczian coordinates (x, y, z) by the spherical polar coordinates (r, 0, (p) in the expressions for the potentials 
obtained from the Lienard-Weichert potential. We have mentioned earlier that the potentials so obtained 
do not satisfy the wave equations in the spherical polar coordinates. 

An interesting aspect of the potential of an oscillating point charge obtained from the distributions (127...) 
is the following. The (i) terms in the higher order normal modes fi{r,r') contributes {■^) radiation terms to 
the electric potential with angular frequency U! where u) is the frequency of the oscillating point charge. Sim- 
ilar discussions are also valid for the corresponding magnetic vector potential. In the dipole approximation 
the amplitudes of these terms are of higher order smallness in the parameter (ks). 

We now construct the solution of the homogeneous wave equation with r < r'. Where (r) is the radial 
coordinate of the field point and (r') is the radial coordinate of the source point. Considering the solutions 
of the static Laplace's equation and the solutions of the wave equation for r > r' we assume that the radial 
solutions of the wave equation for r < r' are given by expressions of the following form: 



Xi = oe{r,l)r''exp{ikr) (141) 

If we substitute the above expression in the radial part of the wave equation we have the following equation 
for the function a{r,l): 

ra" + [2{l + 1) + 2ikr]a' + 2ik{l + l)a = (142) 

It is easy to solve the above equation and the complete potential for r < r' is given by the real or imaginary 
part of the following expression depending on the nature of the source: 



I 3 ^ 

Here the summations over {I) and (j) are running from zero to infinity and if we solve the indicial equation 
[46] for arbitrary values I together with the required regularity conditions at the origin we have the following 
recurrence relation, 



(^ + i + i) 

^(j + l)(j + 2Z + 2)^ 



= -(2z/c)[ ,, , \ ';,'/, ,J a,- (144) 



we chose N to be given by the following expression: 



N = (145) 
47reo 



The constants ao{l) are determined from the source similar to the situation of r > r' . For a source 
oscillating according to a cosine function of time the above solution reduce to the corresponding static solution 
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of the Laplace's equation in the static hmit of static source and c oo if we choose ao{l) = J pQ{r')dv' and 

put the r' dependent terms into the integral Here po{r') is the amphtude of the oscillating source. 

We now demonstrate that in terms of an arbitrary curvilinear coordinate system it is in general not possible 
to find a gauge so that the potentials satisfy both the wave equations and the Lorcntz gauge condition. Let 
be a set of potentials which satisfy the wave equation but d^Af^ ^ 0. If there exists a function A such 
that A^ + d^X satisfy both the Lorentz gauge condition and the wave equations then we have the following 
two conditions on A respectively: 

V^A-10 = -a,A'' (146) 

and 



^*[V^A-^i^]=0 (147) 

The above two conditions are obtained from the Lorentz gauge condition and from the fact that the electric 
potential satisfies the wave equations. The above two conditions are in general contradictory. This happens 
to be the situation with the fields of an oscillating electric dipolc. 

In this section we have only considered the radiation field from oscillating charges. A more complete de- 
scription of radiation problems can be possible through an appropriate quantum theory with the distributions 
determined from a quantum theory. The discussions in this article may also indicate that in all physically 
allowed space-time there exist a preferred coordinate system where the Poisson's equation is satisfied for a 
point charge and also the electromagnetic potentials and the electromagnetic fields are conjugate harmonics 
to each other in four dimension in the Lorentz gauge. This appears to be the situation in the flat space-time 
with the Cartezian coordinate system. 



XXI. MISCELLANEOUS 



One can observe Newton's rings with plane-surface of the convex lens on the glass plate. 

We now discuss an elementary method to calculate he surface integral of a vector field. Let us consider 
the plane ^ + | + | = lin the positive octant. We consider curves on the plane parallel to the X — Z plane. 
The curves are given by the following equation: 

- + - = !-! (148) 
a c 

where y acts as a parameter. An infinitesimal line element on these curves arc given by du = {dx)i — 
{^){dx)k. Similarly we can consider curves on the plane parallel to the Y — Z plane. The curves are given 
by the following equation: 




where x acts as a parameter. An infinitesimal line element on these curves are given by dv = {dy)j — 
{f){dy)k. An area element on the plane is given by du x dv and is expressed in terms of dxdy. We can 
thereafter do the surface integral. One can generalize this method properly. These discussions may also be 
relevant to find out whether a scalar field has an equipotential surface. 

In the following section we will describe the force configuration acting on a uniform cube with weight 'W 
resting on a horizontal plane with friction. The coefficient of static friction is '/i' and the sides of the cube 
are of length '2a'. We consider a coordinate system on the surface with the origin, O, directly below the 
center of the cube. Let a force 'P' act on the block which is not large enough to move the block. The force 
acts along the X-axis. Since the block is not moving the friction between the surfaces should develop an 
equal and opposite force ' — P' acting on the cube at the interface balancing 'P'. However if we consider 
the torque balance w.r.t the origin then the reaction force at the interface will not be uniform. The reaction 
force per unit area 'R' becomes a function of x such that: 
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Pa = 2a f xR{x)dx (150) 

J —a 

The weight is given by a surface integral of R{x) and we have the following two equations: 

P = 2 f xR{x)dx (151) 

J —a 

and 

W = 2a f R{x)dx (152) 

J —a 

It is expected that R{x) will be of the form: R + r{x, P), where R is constant so that the surface integral 
of i? is and r{x, P) is an antisymmetric function w.r.t x such that the equation just above remains valid. 
When P is large enough so that the block is on the verge of motion we have the following relation: 

Pth = 2iia f Rth{x)dx (153) 

J —a 

The first equation then leads to the following condition on r(a;, P): 

nW = 2f xr{x,Pth)dx (154) 

J —a 

which may be useful to determine r{x, Pth) if we consider linear relationships. 

XXII. END 

Project Euclede (Section:ll) is complete. 

Project Galileo-Maxwell (Supplement:II,III) is complete. 
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